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INTRODUCTION 


As  compared  to  the  work  oo  axial  machinery,  very  Little  nos tc  experiments, 
work  has  beet.  done  oo  mixed  and  racial  *ype  compressors,  particularly  for 
nun  speed,  high  pressure  ratio  machines .  This  is  due.  primarily,  to  the 
fact  that  axial  compressors  have  exploited  greater  pro nxse  for  application 
to  aircraft  propulsion.  Other  factors  are:  (1)  the  relative  compactness  of 
the  radial  machine  in  which  the  internal  aerodynamic  processes  are  more 
difficult  to  dissociate  than  in  the  multistage  axiai  machines.  (2)  the  aerodynam. 
complexity  of  the  radial  machine  and  the  greater  numoer  of  interacting 
variables  which  enormously  complicate  systematic  investigations  simuar 
to  those  achieved  on  stationary  cascades.  (31  the  nature  of  the  internal  flow 
in  which  centrifugal  and  Coriolis  effects  play  major  roles  lioth  on  the  main 
flow  and  the  boundary  layer,  thus  rendering  theoretical  compctai*oos  difficult 
and  stationary  testing  useless. 

While  the  axial  compressor  remains  unchallenged  for  tar  bo  mac  nine 
applications  recutring  medium  to  large  through  flows  with  mini  mam  over- all 
diameter,  the  mixed  and  radial  compressors  are  first  choices  for  a  great 
many  other  specific  applications.  This  choice  may  oe  dictated  by  oosKie ra¬ 
tions  of  specific  speed,  limited  axial  length,  ruggedness  simplicity  or  the 
need  of  operating  in  a  dusty  environment .  iocfi  specific  applications  are 
innumerable,  embracing  the  fields  of  aircrait  accessories,  auxniary  power 
plants,  small  aircraft  engines,  ground  equipment,  etc.  In  onef.  tar a> machines 
of  email  to  medium  power  outputs.  -Tie re  the  over- all  diameter  is  not  a 
limiting  factor,  fall  into  this  category.  The  need  for  applied  research  aimed 
at  providing  an  increased  understanding  of  tie  flow  phenomena  in  rr.ixeo  and 
radial  compressors  is  obvious,  if  smaller.  Lighter,  and  mare  dependable 
equipment  is  desired.  Such  research  must  be  a^med  at  providing  mining 
design  information  and  improving  the  performance  from  the  Siam  pc  mis  i 
efficiency,  width  of  the  operating  range,  attainable  pressure  levels  and 
through  flow  capacities.  The  need  for  this  wore  is  made  more  imperative 
by  tr>e  inherent  complexity  of  the  flow  phenomena.  In  addition,  progress 
made  and  the  techniques  developed  in  connection  with  axial  machines  are 
available  and  will  be  of  great  value  reward  improving  our  understanding  of 
the  radial  machine. 

To  this  date,  test  work  on  mined  and  radial  flow  machines  to  investigate 
experimentally  the  effect  of  changes  of  geometry  has  oeen  done  mamiy  by 
NACA .  The  information,  however,  was  not  obtained  fer  the  specific  purpose 
of  establishing  design  criteria.  The  original  intention  was  either  to  investigate 
simple  and  inexpensive  modifications  m  the  hope  of  achieving  some  performance 
improvement,  or  to  alter  a  region  of  the  mpelier  where  some  design  uncertain¬ 
ty  existed  unti*  xe  optimum  shape  was  arrived  at  experimentally. 
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In  references  (36)  and  (37),  the  amount  of  deceleration  in  the  impeller 
was  reduced  successively  In  two  ways,  first,  by  thickening  the  blades,  and 
second,  by  narrowing  the  exit  width.  While  measurements  taken  inside  the 
impeller  show  an  increase  in  efficiency,  this  improvement  is  net  reflected 
in  the  over-all  performance  of  thJ*  compressor  because  of  extraneous  effects 
such  as  changes  in  the  thickness  of  the  blade  wane  and  changes  in  leakage. 
These  results,  coupled  with  the  fact  that  the  lest  imoeller  was  run  ai:  very 
low  tip  speeds,  detract  considerably  from  the  usefulness  of  the  basic  infor¬ 
mation  for  design  purposes. 

In  references  (2)  to  (7),  experimental  study  was  made  of  two  families 
of  impellers .  The  two  families  differ  m  that  the  flow  area  is  -ontrolltd 
in  one  case  by  the  meridional  shape  and  in  the  other  by  the  olade  thickness 
variation.  Similar  blade  loadings  were  aimed  at  but  were  not  achieved. 

Within  each  family,  variation  of  the  tip  exit  width  was  introduced  to  change 
the  amount  erf  eieceleration  in  the  impeller.  The  test  results  clearly  brought 
out  the  superiority  oi'  the  thin-bladed  impeller,  although  one  can  only  gue*s 
at  the  causes.  Difficulty  also  exists  in  evaluating  the  effect  of  changing  tip 
width,  which  is  not  clearly  brought  out  by  test,  probably  doe  to  simultaneous 
cianges  in  the  diffuser  performance. 

The  effect  of  blade  loading  was  investigated  in  references  (2i»s  (22),  and 
other  reports  yet  to  be  published .  Preliminary  information  available  at  this 
date  seems  to  indicate  that  better  performance  will  be  achieved  with  a  lightly 
loaded  inducer. 

It  is  seen  that  there  is  only  meager  systematic  investigation  of  the 
basic  parameters  in  the  field  of  radial  or  mixed -flow'  compressors.  Without 
this  information,  basic  design  criteria  cannot  be  set  and  fundamental  im¬ 
provements  in  compressor  performance  are  left  to  chance.  It  is  the  prime 
objective  of  the  present  program  to  isolate  a  few  more  basic  and  significant 
variables  in  order  to  investigate  them  sovb”,  to  evaluate  their  effects  upon 
the  impeller  performance,  and  ultimately  to  set  up  basic  criteria  for  the 
design  of  high  performance  compressors. 


xvi 

CONFIDENTIAL 


WADC  TR53-589 


CONFIDENTIAL 


I  STUDY  OF  VARIABLES  AND  DESIGN  CRITERIA 

The  first  objective  of  the  present  research  program  Is  to  select  among 
ail  the  variables  erMrLTg  Sie  d£t*gc  of  2n  .mpeller  those  which  first, 
are  universal  enough  to  be  used  for  any  impeUer  design,  and  second,  are 
of  greatest  significance  ui  their  effect  upon  the  impeller  performance,  rte 
second  objective  of  this  program  ls  to  investigate  theoretically  and  experi¬ 
mentally  these  effects  to  establish  basic  design  criteria. 

A.  Discussion  of  Design  Variables 

Within  the  framework  of  externally  .mposed  impeller  design  features, 
based  on  non-aerodynamic  require  ^ects  cr  on  aerodynamic  considerations 
the  effects  of  which  are  thought  to  be  well  -nderstacd.  there  remains  a 
considerable  amount  of  design  freedom,  particularly  in  the  internal  shape 
erf  the  impeller.  Those  aerodynamic  features  which  are  be^ered  to  be 
understood  become  part  of  the  regular  design  procedure.  They  are  folly 
described  in  Section  II.  General  Impeller  Design  Procedure.  The  remaining 
aerodynamic  features  required  to  complete  the  design  are.'  for  the  most 
part,  established  from  experience  and  judgement  As  will  be  explained 
they  consist  mainly  of  the  choice  of  the  meridional  contour,  the  Dumber, 
shape,  and  thickness  of  the  blades,  the  distribution  and  amount  of  incidence, 
and  the  length  of  the  flow  passages.  It  is  among  these  variables  that  selectman 
of  the  most  significant  quantities  must  be  made  in  the  present  study.  Only 
a  few  of  these  variables  can  be  considered  secaraieiv.  if  the  total  number  of 
variations  is  to  be  kept  within  reason.  These  few  will  be  called  primary 
variables.  The  question  then  arises  as  to  what  constant  value  to  assign  to 
the  so-called  secondary  variables.  One  approach  is  that  the  secondary 
variables  be  optimized  within  the  design.  Each  design  would  then,  m  itself, 
be  an  optimum  configuration  for  the  particular  variation  coosen  .  The  problem 
ultimately  would  then  be  reduced  to  the  finding  of  the  best  erf  several  optimized 
impellers  or,  figuratively  spraxing.  finding  the  highest  peak  among  several 
peaked  curves.  However,  it  is  quite  apparent  that,  if  such  an  approach  were 
rigidly  followed,  a  good  impeller  would  undoubtedly  be  produced  but  it  would 
be  virtually  impossible  to  isolate  the  separate  effects  of  the  primary  variables 
because  of  variation  in  the  secondary  variables  Another  approach  would  be 
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to  keep  the  secondary  variables  unchanged.  This  has  the  disadvantage  of  not 
being  able  to  produce  an  optimum  impeller,  because  the  ultimate  optimum 
combination  of  primary  variables  would  not  be  achieved  with  the  optimum 
combination  of  secondary  variables  (which  have  remained  unchanged/.  It 
appears  then,  that  a  compromise  approach  becomes  necessary.  Optimization 
of  each  primary  variable  should  be  carried  out  only  so  far  as  variation  of 
the  secondary  variables  is  not  required.  Each  primary  design  varlatior  must 
be  so  made  tnat  all  secondary  variables  can  remain  unchanged  and  yet  make 
up  a  near  optimum  design.  If  this  compromise  is  practical,  it  will  hr.ve 
the  advantage  of  isolating  the  significant  variables  and  sull  produce  near 
optimum  impeller. 

B.  Selection  of  Design  Variables 

It  was  quite  apparent  from  the  start  that  the  choice  erf  the  design  variation 
should  be  limited  to  variables  internal  to  the  impeller  rather  than  tc  variable 
relationships  between  the  impeller  and  its  external  surroundings,  r>uch  as 
mean  line  velocity  triangles,  impeller  cone  angle,  inlet  geometry,  and  the  like 
This  will  be  done  consistently  in  this;  work,  except  as  is  otherwise  specified 
in  the  original  research  contract. 

1 .  Over-all  Blade  Loading 

One  of  the  most  fundamental  quantities  entering  the  design  of  the 
compressor  in  general  is  the  over -all  blace  loading.  This  criteria  is, 
at  the  present  time,  missing  in  centrifugal  impellers  and  is  left  to  the 
experience  of  the  designer.  What  is  needed  is  the  counterpart  of  such 
axial  compressor  design  quantities  as  the  difiusion  factor  or  the  product 
of  solidity  times  lift  coefficient.  This  need  is  particularly  important 
in  view  of  the  fact  that  the  radial  nature  of  the  flow  may  permit  the 
achievement  of  quite  high  loading,  if  advantage  is  takeu  of  potential 
ooundary  layer  control  by  centrifugal  action. 
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Achievement  of  variable  over -all  blade  loading  can  be  done  either 
by  changing  the  number  of  blades  or  by  changing  the  impeller  flow  length, 
each  affecting  directly  the  blade  solidity.  Changing  of  the  number  of 
blades  is  not  altogether  satisfactory  m  view  of  the  fact  that  in  high 
specific  speed  machines  an  upper  limit  to  the  number  of  blades  is  soon 
reached  because  of  restrictions  imposed  by  cutter  size  requirements 
in  machining  the  blades  in  the  hub  legion.  A  more  flexible  solution  is 
to  vary  the  over -all  blade  loading  by  combining  changes  in  the  number 
of  blades  together  with  changes  in  blade  chord,  i.e. ,  changes  in  impeller 
length.  It  i3  difficult  to  predict  which  of  the  two  effects  is  best  from 
the  standpoint  of  efficiency.  Wall  friction  is  about  the  same  for  toth. 

The  extent  of  the  secondary  flow  in  both  cases  is  cortr oiled  by  the 
pressure  gradients  in  the  boundary  layer,  the  magnitude  and  directioa  of 
which  is  affected  by  many  considerations  besides  passage  length  and 
number  of  blades.  In  favor  of  the  impeller  with  many  blades,  and  a 
short  passage  is  the  fact  that  secondary  flow  accumulations  are  quietly 
discharged  through  the  exit,  but  this  action  is  somewhat  offset  by  the 
fact  that  greater  design  flexibility  exists  with  longer  passage  sad  fewer 
blades.  Again,  in  favor  of  the  former,  is  the  lighter  weight  cl  the  impeller. 

It  appears  then  that  an  impeller  'with  iightiy  loaded  blades  will  have 
a  general  appearance  of  being  longer  and  will  offer  greater  design 
flexibility  than  its  stubbier  counterpart  with  heavily  loaded  blades .  It 
must  here  be  brougnt  out  that  the  success  ot  failure  cf  either  type  will 
have  profound  implication  upon  the  form  of  the  best  design  philosophy . 

The  former  approach  ^lightly  loaded  blades  and  long  passages)  is  tased 
upon  a  design  philosophy  in  which  large  frictional  effects  are  tolerated 
but  where  boundary  layer  build-ups  are  minimized.  The  latter  (heavily 
looded  blades  and  short  passages)  iS  based  on  the  idea  that  some  flow 
separation  is  anticipated  but.  because  of  the  shortness  of  the  impeller 
and  the  shape  of  the  meridional  section,  stagnating  fluid  will  be  dischar^ 
before  it  can  build  up  in  one  area  and  d'srupt  the  over -ail  flow  pattern. 
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2 .  Distribution  of  Blade  Loading 

Just  as  important  as  the  amount  of  over-all  loading  t«  the  distri¬ 
bution  of  that  loading  between  inducer  end  and  wheel  outlet  end.  A 
heavily  loaded  inducer  ie  attractive  in  that  considerable  work  is  done 
where  the  boundary  layer  is  new  and  little  developed.  However,  it  is 
felt  that  because  of  the  high  sensitivity  of  the  flow  in  a  transonic  inducer 
the  above  achievement  is  difficult.  A  heavily  loaded  wheel  portion 
suffers  from  accumulated  boundary  layer  and  from  developed  secondary 
flow,  but  may  be  success! 'll  if  full  advantage  is  taken  of  the  centrifugation 
of  the  boundary  layer . 

Variation  in  blade  loading  distribution  along  the  impeller  can  be 
achieved  by  changing  the  meridional  shape  or  by  adjusting  the  blade 
thickness  or  blade  camber .  Of  these,  blade  thickness  changes  are  uot 
satisfactory,  since  an  appreciable  fraction  of  the  flow  relative  to  the 
blade  is  supersonic .  In  order  to  minimise  blade  losses,  it  hs  necessary 
to  keep  the  blades  as  thin  as  feasible.  Changing  the  meridional  contour 
to  achieve  changes  in  loading  is  subject  to  two  difficulties.  First,  any 
desired  loading  cannot  be  so  attained,  d  blade  camber,  inlet,  and  exit 
conditions  are  fixed.  Such  contour  changes  are  much  mnri  suitable  as 
a  tool  for  making  refinements  of  a  local  nature  rather  than  basic  changes 
for  the  whole  impeller.  Second,  changes  in  meridional  contour  almost 
always  entail  changes  both  at  shroud  and  hub.  This  causes  loading 
readjustments  between  hub  anc*  shroud,  changes  in  charnel  aspect  ratio, 
etc . ,  with  the  net  result  that  the  effect  of  loading  distribution  between 
inlet,  and  exit  cannot  be  isolated.  A  more  direct  approach  then  is  to 
achieve  the  desired  effect  by  changing  the  shape  of  the  blade  camber  line 
only .  Loading  is  very  sensitive  to  this  latter  variable  and  no  other 
extraneous  effect  is  introduced. 
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3.  Impeller  Exit  Cone  Angle 

Investigation  of  the  effect  o i  varying  the  exit  cone  angle  was 
specifically  requested  in  the  original  problem  statement  for  this 
program .  Even  though  exit  coae  angle  is  an  external  impeller 
variable,  ls  main  effect  is  that  of  varying  the  impeller  length. 

As  such,  it  may  be  considered  as  an  additional  variation  under 
part  (1)  above,  in  addition,  it  is  anticipated  mat  a  redistribution  of 
blade  loading  between  shroud  and  hut)  will  take  place.  Ihe  hub 
loading  in  the  inducer  will  t>e  even  lower  in  the  case  of  the  30-degree 
impeller  than  for  the  60-degree  txit  case,  first,  because  of  the 
lesser  fluid  deceleration  there  and  second,  because  turning  of  the 
olade  at  the  hub  (with  radial  blade  elements!  will  be  less  severe  than 
for  a  60-degree  case  with  the  same  blade  shape  ai  the  sbroun 

At  the  present  ome,  it  would  be  extremei>  hazardous  to  attempt 
to  make  predictions  of  the  efficiency  potentialities  of  the  30-degree 
impeller  as  compared  to  the  60-degree  type.  It  may  be  said,  however, 
that  a  good  part  of  the  answer  may  be  Known  when  the  results  of 
part  (1)  above  are  available. 

it  must  be  Kept  in  mind  that  an  impeller  with  a  20 -degree  cooe 
angle  will  require  a  different  diffuser.  As  such,  the  potentialities 
of  the  compressor  must  not  be  confused  with  the  potentialities  of  the 
impeller.  In  addition,  the  impeller  with  a  low  exit  cone  angle  will 
almost  necessarily  be  a  heavier  piece  of  machinery. 
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n  GENERAL  IMPELLER  DESIGN  PROCEDURE 

Contrary  to  the  aerodynamic  design  information  available  for  axial-flow 
compressors,  no  generally  accepted  procedure,  backed  by  extensive  fuau 
detailed  teat  information,  exists  for  radial  or  mixed-flow  machines.  Partly 
because  of  the  complex  nature  of  the  flow  process  inside  su;h  impellers 
and  because  of  the  large  number  of  possible  geometrical  parameters,  no 
known  systematic  testing  U  impellers  or  impeller  components  has  been 
undertaken.  Whatever  information  is  available  on  such  things  as  details 
of  the  internal  flow  patterns  or  the  effect  of  some  particular  geometrical 
variation  la  restricted  either  to  liquid  pumps  cr  to  low  speed  centrifugal 
compressors,  (ref.  35.37.)  Design  of  the  higher  speed  compressors  is 
based  mainly  on  theoretical  considerations  combined  with  whatever 
Information  can  La  borrowed  from  the  wor1:  done  will,  axial  compressors  and 
with  centrifugal  compressors  handling  incompressible  or  nearly  incompressible 
fluids.  In  short,  experimentally  determined  design  criteria  are  net  «*.J,i»bie 
to  supplement  theory  in  those  areas  where  the  theory  is  Insufficiently  deveto,  ed 
to  supply  the  information 

Among  the  many  such  design  criteria  may  be  listed:  (1)  the  best  loading 
distribution  between  inducer  and  radial  portion  of  impeller,  (L'  a  means  of 
predicting  and  controlling  the  behavior  of  the  bounds  y  layer,  (3)  an  opH*num 
combination  of  blade  turning  and  meridional  deceleration,  (4)  the  effects  o> 
the  loading  distribution  between  hub  ami  shroud,  (5)  which  loading  (bub,  akrosa 
or  50  per  cent  streamline)  must  be  controlled  most  carefully  in  that  it  bears 
greatest  influence  on  the  over-all  performance  of  the  compressor,  (6)  a 
quantitative  understanding  of  the  sources  of  losses,  their  magnitude  aqd 
di  sir i  notion,  (7)  the  effect  of  impeller  geometry  on  the  shape  of  the  exit 
velocity  profile,  etc.  In  the  absence  of  such  general  information,  otbek 
criteria  must  be  devised.  Some  are  based  on  a  meager  understanding  of  the 
internal  flow  behavior  from  theoretical  considerationc  ad  others  are  bused 
on  judgement  and  experience. 
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It  is  generally  felt  that  a  good  design  must  not  only  be  such  that  all 
calculated  flow  quantities  vary  in  a  manner  acceptable  to  theory  but  al*o  that 
the  design  as  a  whole  must  be  appealing  to  the  experienced  eye.  This 
restriction  rules  out  oddities  in  design  such  as  startling  wall  contour  shapes, 
wanr  j  flow  passages  etc .  It  is  also  felt  that  the  greater  care  in  design  must 
be  given  to  the  shroud  rather  than  to  the  hub .  This  is  supported  not  only  by  the  :  act 
that  the  larger  fraction  of  the  total  flow  is  subjected  to  the  shroud  conditions, 
but  also  because  relative  velocities  there  are  higher.  Low  14ach  number  flews 
Lave  the  ability  to  adjust  themselves  comparatively  easily  to  geometrical 
variations,  whereas  h?gfa  fcurh  number  flows  are  extremely  sensitive  to  such 
changes.  In  addition,  losses  incurred  in  high  velocity  flow  are  usually  large. 

At  thi£  time,  indications  are  that  the  variation  of  suction  and  pressure 
side  relative  velocities  is  a  good  index  of  the  quality  of  a  design.  Difficulties 
arise  only  when  the  criteria  of  rate  of  deceleration  and  amount  and  distribution 
cf  loading  are  needed.  In  the  absence  of  such  specific  data,  we  must  satisfy 
ourselves  with  the  following  rules ; 

(1)  All  flow  passage  areas  must  exhibit  smooth  and  gradual 
variations . 

(2)  All  velocities  and  pressures  along  streamlines  must  change 
in  a  smooth  continuous  way,  with  avoidance  of  high  localised 
gradients. 

(3)  Avoidance  of  ell  unnecessary  accelerations  or  decelerations . 

If  an  over-all  deceleration  process  is  desired,  for  example, 
the  intermediate  process  must  be  decelerating  everywhere, 
except  for  unavoidable  exreptions  where  accelerations  must 
be  kept  at  a  minimum . 

(4)  Relative  velocity*  must  be  kept  comfortably  positive  everywhere 
in  order  to  avoid  dead  fluid  areas  or  back  flow. 

( 5)  Impeller  blade  loading  must  be  uniformly  distributed,  with  the 
loading  limit  on  the  inducer  not  greatly  exceeding  that  attainable 
in  a  corresponding  axial  cascade 
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(6)  Large  decelerations  in  regions  of  anticipated  thick  bounds*  3 
layers  must  be  avoided.  Reference  (1)  provides  information 
on  boundary  layer  growth  against  adverse  pressure  gradients. 

(7)  A  certain  amount  of  control  of  secondary  flow  can  be  had  by 
designing  In  such  a  way  that  the  radial  pressure  gradients 

in  the  main  flow  will  tend  to  neutralise  the  centrifugal  pressure 
gradieuts  set  up  on  the  blades  in  the  boundary  layer . 

In  transonic  flow,  it  is  fell  that  most  of  the  difficulty  experienced  is  due 
to  our  inability  to  describe  the  flow  with  sufficient  accuracy .  Transonic 
flow  is  extremely  rensitive  to  area  changes.  Design  procedures  which  predict 
the  flow  passage  areas  with  sufficient  accuracy  for  low  subsonic  velocities 
become  inadequate  at  high  Mach  numbers  because  approximations  are  no  longer 
tolerable  without  great  sacrifices  in  efficiency.  With  the  development  of  more 
exact  design  procedures  it  is  reasonable  to  expect  that  the  efficiency  of  high 
flow,  high  pressure  ratio  machines  can  be  made  as  high  as  their  more  con¬ 
servative  relatives. 

A.  One- Dimensional  Aspects  cf  Design  and  Determination  of  External 
Geometry  of  Impeller 

It  is  intended  in  this  section  to  present  the  procedure  by  which  the  basic 
external  dimensions  of  an  impeller  are  established.  While  the  procedure  is 
based  on  concepts  which  are  generally  accepted,  a  good  many  different  approaches 
exist .  The  present  one  (considered  here)  has  the  advantage  of  reducing  and 
grouping  several  steps  of  the  procedure  into  a  simpler  form  which  is  applicable 
to  any  radial  or  mixed -Dow  impeller . 

Hie  usual  requirements  of  a  compressor  are  that  it  be  capable  of  handily 
a  specified  mass  rate  of  flow  against  a  specified  pressure  ratio.  This  informa¬ 
tion  together  with  some  knowledge  of  the  particular  application  of  the  unit  estab-ish 
almost  completely  the  one -dimensional  over -all  features  of  the  compressor .  < 

A  small  amount  of  freedom,  however,  is  left  to  the  designer.  It  is  one  of  the 
purposes  of  the  present  research  program  to  single  out  these  degrees  of 
freedom  with  the  intention  of  establishing  criteria  leading  to  optimum  designs. 
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It  is  necessary  to  differentiate  between  impeller  external  and  internal 
geometry.  In  the  external  geometry  there  are  included  all  quantities  and 
dimensions  which  are  independent  of  the  details  of  the  flow  pattern  inside 
thr  impeller  passages.  The  external  geometry  can,  in  general,  be 
established  by  one- dimensional  study  of  the  inlet  and  exit  conditions.  It 
involves  the  following  quantities: 

Impeller  blade  angle  at  the  exit 

Impeller  tip  speed 

Impeller  mean  exit  diameter 

Inducer  tip  diameter 

Inducer  huo  diameter 

Impeller  exit  width 

Impeller  (approximate/  over-all  length 
Impeller  exit  cone  angle 

On  the  other  hand,  the  internal  geometry  of  the  impeller  can  only  be 
arrived  at  by  detailed  two-  and  three-dimensional  analysis  of  the  inter^c.1 
flow.  The  Internal  geometry  is  completely  defined  by  the  following  quantities. 

Exact  shape  of  hub  and  shroud  contour 

Blade  camber  distribution 

Humber  of  blades 

Blade  thickness  distribution . 

The  details  of  the  analysis  of  the  internal  flow  will  be  the  subject  of  the 
next  chapter . 
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As  a  rule,  it  is  not  possible  to  solve  directly  for  each  external  quantity 
In  succession.  The  procedure  is  rather  oue  of  successive  approximations  to 
be  repeated  until  all  conditions  are  met  and  all  the  various  physical  QsiatUiM 
are  consistent  throughout. 

Two  of  the  above  listed  internal  quantities  are,  at  the  present  time, 
without  benefit  of  specific  design  criteria  and  are  left  to  the  designers 
experience:  they  are  Impeller  length,  and  Impeller  exit  cone  angle.  As 
previously  mentioned,  these  quantities  have  been  selected  as  two  of  the 
primary  variables  to  be  studied  In  the  present  program.  It  is  believed  that 
all  other  external  physical  quantities  are  fixed  by  specific  design  criteria. 

A  discussion  of  these  criteria  will  now  be  undertaken;  taken  together  they 
make  up  the  design  procedure  followed  by  AlResearch  for  tbs  Impeller 
external  features . 

(1)  Impeller  blade  angle  at  exit. 

In  reference  (1)  a  generalized  study  was  undertaken  to  optimize  the 
geometry  of  radial  impellers.  The  results  are  applicable  to  impeller  a  with 
slightly  mixed-flow  as  well.  The  study  Indicates  that,  for  the  type  of 
impellers  here  considered,  the  effect  of  blade  exit  angle  on  over-all 
compressor  efficiency  is  slight  provided  the  blades  are  nearly  radial. 

In  addition,  the  Impeller  with  nearly  radial  blades  was  found  to  bo  optimum 
in  efficiency  for  the  range  of  pressure  ratios,  here  considered. 

A  further  consideration  entering  Into  the  selection  of  blade  exit  angle  Is 
that  of  stress.  High  Up  speed  dictates  the  choice  of  blades  consisting  of 
radial  elements.  In  strongly  mixed-flow  impellers  backward -bent  blades 
can  be  had  with  radial  elements  but  for  Impellers  which  are  only  slightly 
mixed-flow  (nearly  radial)  the  amount  of  backward  curvature  possible  is 
diminiehingly  small. 

(2)  ImpeUer  Up  speed. 
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The  selection  erf  the  impeller  tip  speed  depends  on  the  choice  o'  the  mean 
inlet  and  exit  velocity  triangles  These  triangles. in  turn,  must  satisfy  the 
requirement  oi  impeller  work  input.  While  the  amount  orf  work  input  is 
primarily  dictated  by  the  required  pressure  ratio,  it  also  involves  considera¬ 
tion  of  the  extent  and  distribution  of  the  losses  in  the  compressor.  The 
assignment  of  numerical  values  to  these  losses  must  obviously  be  based  on 
experience  and  available  test  results.  This  circumstance  holds  equally  true 
for  several  other  factors  required  to  complete  the  velocity  triangles.  They 
are- 


(1)  Amount  of  impeller  slip  ( ref.  1,  8,  9,  10.  11) 

(2)  The  radial  velocity  gradients  anticipated  along  the 
inducer  leading  edge. 

(3)  Maximum  permissible  relative  flow  deceleration  in  the 
impeller,  expressed  either  in  the  f'jrm  of  a  chifusian  factor, 
or  more  commomy,  as  a  deceleration  ratio,  (ref.  10.  11) 

Figure  1  summarizes  some  of  the  results  of  ref  (1).  As  might  be  antici¬ 
pated,  the  impeller  efficiency  is  strongly  affected  by  the  amount  of  wheel 
reaction,  i.e. ,  high  values  erf  the  deceleration  ratio  W^/W,  produce  high 
impeller  efficiencies.  Less  obvious,  however,  are  the  indications  that  the 
compressor  efficiency  also  improves  for  decreasing  wheel  reaction.  This 
trend  is,  in  general,  also  supported  by  a  study  of  a  large  number  of  compressors 
in  actual  operation.  It  is  odvious.  moreover,  that  a  limit  of  efficiency  must 
soon  be  reached,  since  for  zero  reaction,  the  ccmpi essor  efficiency  will 
be  below  that  orf  an  overloaded  diffuser.  Furthermore,  it  is  clear  that  the 
over-all  size  orf  the  unit  must  aiso  increase,  if  a  reasonable  diffusion  rate 
is  to  be  maintained.  The  optimum  amount  of  impeller  reaction  is  a  design 
variable  of  great  importance  and  could  well  form  the  subject  orf  an  independent, 
extensive  research  program.  Information  on  the  effects  of  specific  speed 
and  compressibility  for  an  optimum  degree  of  reaction  are  absolutely  needed. 

In  the  absence  of  such  data,  one  usually  adopts  a  compromise  value  of 
impeller  reaction  at  from  55  to  60  per  cent,  (a  value  above  50  per  cent  being 
selected  in  view  orf  the  fact  that  boundary  layer  conditions  are  felt  to  be  more 
favorable  in  the  impeller  than  in  the  diffuser).  In  practice  most  orf  these 
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factors  require  re-cvaluation  as  the  design  progresses  and  more  of  the 
impeller  geometry  becomes  known.  These  modifications  can  usually  best 
be  accommodated  by  slightly  changing  the  exit  blade  angle. 

(3)  Impeller  Mean  Exit  Diameter. 

Since  the  compressor  flow  and  total  head  are  prescribed,  selection  of 
the  impeller  mean  exit  diameter  establishes  the  impelier  specific  speed 
(equation  A  -  5),  which,  in  turn,  very  eiosexy  sets  most  of  the  remaining 
impeller  external  dimensions.  Selection  of  the  impeller  diameter  is  often 
made  by  considerations  other  than  purely  aerodynamic  ones:  (1)  since  the 
selection  determires  the  physical  speed,  it  must  be  so  made  as  to  provide 
favorable  matching  between  turbine  and  compressor .  In  a  complete  unit 
design,  therefore,  the  selection  of  impeller  diameter  is  in  variably  tied 
to  the  turbine  design.  Obviously,  for  the  purpose  of  the  present  program, 
a  somewhat  greater  degree  ol  freedom  is  permissible,  since  the  turbine 
matching  requirement  is  reduced  to  the  requirement  of  tailoring  the  com¬ 
pressor  to  fit  an  aready  existing  convenient  driving  unit.  (2)  A  further 
consideration  is  that  of  physical  size  and  weight. 

Small  envelope  dimensions  together  with  lightweight  impose  conditions 
of  high  rpm  and  high  absolute  and  relative  flow  Mach  numbers,  i.e.,  a 
critical  design.  These  requirements  token  together  completely  define  the 
impeller  tip  diameter . 

(4)  Inducer  Tip  Diameter . 

A  convenient  expression  for  the  ratio  of  inducer  tip  diameter  to  mean 
exit  tip  diameter  can  be  deduced  by  making  use  of  the  concept  of  specific 
speed  and  by  designing  for  a  minimum  inducer  tip  relative  Mach  number . 
Combination  of  expressions  A-l  and  A -5  for  specific  speed  yields  the 
following  convenient  expression 


_j _ r  _u_  Kf  i  i'*  D  .  Ad~ 

6.69  Wcy^  ton  £s  J  inches  V  ds 
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which  is  derived  in  Appendix  A-l-d.  In  this  relationship,  everything  is 
known,  except  the  value  of  $s-  This  angle  must  be  varied  to  achieve  the 
optimum  condition  of  minimum  inducer  tip  relative  Mach  number.  This 
optimization  was  carried  out  ir.  Appendix  A- 1  -b,  in  which  it  is  shown  that 
the  optimum  angle  is  a  function  of  mean  inlet  Mach  number  only  (see 
Figure  2).  Thus,  for  incompressible  flow,  the*  angle  for  minimum  relative 
velocity  is  54.8  degrees  For  a  mean  inlet  Mach  number  of  0  5.  this 
*ngie  is  58.5  degrees.  Combining  these  conditions  with  equation  H- 1-1 
above  yields  the  relationship  illustrated  on  Figure  3.  From  this  figure, 
the  desired  inducer  diameter  may  be  obtained  irom  the  known  values  of 
compressor  flew,  tip  speed,  and  exit  diameter,  and  from  estimates  of 
inlet  mean  Mach  number,  inlet  velocity  gradient,  and  inducer  hub-tip 
ratio.  The  possible  error  introduced  by  misjudgement  of  these  various 
estimates  Is  small,  first  because  the  inducer  diameter  is  a  weak  function 
of  Lnese  variables  and  second,  because  near  the  optimum  point,  small 
variations  In  inducer  diameter  result  in  insignificant  changes  in  inducer 
tip  relative  Mach  .lumber . 


(5)  Inducer  Bub  Diameter. 

From  considerations  of  unit  size  and  weight,  it  is  obvious  that  small 
inducer  bub  diameters  are  desirable  in  order  to  increase  the  gulp  capacity 
of  the  unit.  Restrictions  are  Introduced  by  the  stress  Limitations  at  the  blade 
roots  and  by  cutter  limitations  m  machining  the  blades.  These  considerations, 
taken  together  for  the  chosen  number  of  blades,  determine  the  minimum 
possible  hub  diameter .  At  this  stage  of  the  design,  the  choice  of  the  number 
of  blades  may  or  may  not  be  final,  however,  the  initial  choice  is  nevertheless 
accurate  enough  for  the  purpose  of  finding  the  hub  diameter . 

(6)  Impeller  Exit  Width. 

The  width  of  the  impt.’ler  at  the  exit  must  be  such  that  continuity  Is 
satisfied  (at  the  design  mass  flow)  for  the  desired  exit  velocity  triangles. 

For  computation  c 4  continuity,  the  state  of  the  fluid  at  the  impeller  outlet 
is  required.  Appendix  A  derives  expressions  for  the  pressure  and  density 
variation  inside  an  impeller  as  a  function  of  the  local  velocity  triangle  and  a 
poly  tropic  exponent  n  .  This  procedure  obviously  assumes  homogeneous 
distribution  of  the  irreversibilities.  Appendix  A  also  gives  a  relationship 
between  the  exponent  n  md  impeller  efficiency. 
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An  increase  in  the  exit  width  by  a  certain  amount,  in  the  fora  of  a 
viscous  clogging  coefficient  *  ,  has  been  used  to  accommodate  the  boimdan 
layer  displacement  thickness.  More  recent  work  has  demonstrated  that 
4)  the  introduction  of  such  an  independent  viscous  (actor  introduces  a 
redundancy  in  the  design  and  cannot  be  justified  theoretically,  (2}  test  data 
indicates  that  the  theoretically  correct  valve  of  «  is  close  to  1.0.  Inasmuch 
as  the  present  design  was  made  to  incorporate  a  variation  of  flow  factor 
as  given  in  Figure  4,  it  is  felt  that  perhaps  slightly  excessive  flow  deceleration 
is  built  into  the  impeller  and  that  the  exit  width  is  too  large.  A  more  elaborate 
treatment  on  the  subject  of  flow  factors  is  developed  in  Appendix  A.  The 
analysis  indicates  that  for  a  reasonably  good  exit  profile,  a  coefficient  of 
the  order  of  0.96  Is  sufficient  for  computing  the  impeller  exit  width. 

(7)  Impeller  exit  cone  angle  and  impeller  over -all  length. 

As  mentioned  earlier,  the  effects  of  these  variables  on  impeller 
performance  are  not  known  and  one  to  be  investigated  as  part  of  this  program. 
The  effects  of  these  variables  on  the  theoretical  quality  c!  the  design  are, 
however,  understood  to  some  extent.  Thus,  for  example,  favorable  shroud 
velocity  distribution  can  be  achieved  only  if  the  design  employs  gentle  shroud 
curvatures.  For  impellers  with  a  high  ratio  of  inducer  to  Up  diameter 
(say  above  0. 5)  gentle  shroud  curvature  can  be  achieved  only  in  an  impeller 
of  the  mixed-flow  type.  Am  approximate  guide  for  shroud  curvature  is  that 
the  mean  shroud  radius  of  curvature  in  the  meridional  view  should  be  equai 
to  about  1.25  times  the  length  of  the  inducer  blade.  This  rule,  in  general, 
results  in  reasonable  velocity  gradients.  It  also  fixes  the  impeller  length. 
Whether  this  approach  can  be  proved  best  experimentally  remains,  however , 
to  be  demonstrated. 

B.  Three-Dimensional  Aspects  of  Design  Procedure  and  Determination  of 
Impeller  Internal  Geometry 

Once  the  exte  rcrd  dimensions  of  the  impeller  have  been  established,  use 
is  made  of  a  combinaUon  of  direct  and  indirect  procedures  for  the  calculation 
of  the  internal  impeller  geometry,  i.e.,  the  meridional  contour  and  the  number 
and  shape  of  the  blades.  The  best  meridional  shape  is  established  indirectly, 
in  that  successive  modifications  are  made  until  all  veiociUes  have  the  desired 
variation  and  all  shapes,  contours,  and  area  variations  have  the  desired 
degree  of  smoothness.  For  each  such  meridional  shape,  the  niade  camber 
and  blade  thickness  are  prescribed  by  a  direct  method  of  solution. 
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AlResearch  is  presently  engaged  in  a  program  o f  confirming  and 
improving  its  compressor  design  procedures.  The  original  procedure  i a 
based  on  the  use  of  hand  computers  and  is  well  adapted  for  design  purposes. 
Necessarily  involved  are  a  number  of  simplifying  assumptions  to  keep  the 
computing  time  within  reason .  An  improvement  of  the  design  procedures 
presently  in  progress  involves  the  carrying  out  on  high  speed  computers 
of  more  exact  analyses  of  the  impeller  internal  flow  for  comparison  with 
the  approximate  solution  and  with  test  data.  If  these  computations  are 
successful,  it  would  be  extremely  desirable  to  apply  this  procedure  to  some 
impellers  of  the  present  program.  Below,  follows  a  description  of  the 
original  design  procedure  and  of  the  more  exact  method  subsequently 
developed . 

1 .  Basic  Design  Procedure 

fa)  Calculation  of  the  Meridional  Flow  in  the  Impeller. 

The  basic  working  approach  of  the  design  procedure  involves 
the  analysis  of  the  flow  through  the  inlet  and  through  the  impeller 
annulus  as  in  vise  id,  irrotationai,  quasi -compressible,  and  quasi- 
three  dimensional  flew .  Corrections  are  mace,  either  concurrently 
or  subsequently,  for  a  more  accurate  description  of  the  real  flow. 

The  corrections  allow  for  viscous  effects,  blade  clogging  effects, 
the  three-dimensional  nature  of  the  flow,  and  for  vorticity  effects 
as  introduced  by  the  nonuniformity  of  the  work  addition  in  the 
impeller . 

The  inviscid ,  irrotationai,  quasi -compressible,  and  quasi - 
three  dimensional  flow  is  taken  as  an  approximate  flow  picture 
showing  all  the  main  properties  of  the  final  flow.  This  flow  description 
is  used  as  a  basic  tool  for  the  study  of  the  characteristics  of 
various  designs.  As  such  it  may  be  defined  as  that  flow  which  has 
the  same  velocity  gradients  along  a  potential  line  from  hub  to 
shroud  as  the  two-dimensional,  incompressible,  inviscid  and 
irrotationai  flow,  but  in  which  satisfaction  of  continuity  is  achieved 
by  using  the  three-dimensional  area  with  the  correct  densities 
corresponding  to  the  calculated  velocities. 
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The  two-aimer/8ional,  incompressible,  invtscid  and 
irrotatlonal  velocity  gradients  are  obtained  by  making  an  analog  field 
plot  of  the  meridional  shape  under  study ,  This  meridional  shape 
includes  both  th  npeller  and  the  inlet  portions.  The  analog  field 
plot  consist  of  an  electrical  analog  cut  from  Teledeltos  conductive 
paper  in  the  trial  meridional  shape.  Hie  equipment  and  the 
procedures  used  are  those  of  ref.  18.  For  the  present  applica¬ 
tion,  a  plot  of  potential  lines  rather  than  streamlines  is  of 
greater  convenience.  The  plot  is  made  by  imposing  a  fixed 
potential  across  two  estimated  conductive  potential  lines,  one 
for  the  upstream  boundary  and  the  other  for  the  downstream 
boundary.  Between  these  boundary  potential  lines,  additional 
equipotential  lines  are  found.  For  incompressible  flow  iind  along 
a  given  potential  line,  their  spacing  Ajp  is  inversely  porportional 
to  the  local  velocity 

Cm  Am 

The  value  of  K  can  be  obtained  by  satisfying  ccntiauity,  as 
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In  this  relationship,  the  quantity  <  is  introduced  as  a  viscous  flow 
factor  and  Zrr  represents  the  blade  blockage  area .  Both  these 
quantities  can  be  prescribed  with  fair  accuracy . 

The  flow  coefficient  *  has  been  computed  experimentally  at 
the  impeller  Inlet  and  exit.  At  the  impeller  inlet,  definition  of 
the  term  *  is  similar  to  that  used  in  a  plain  nozzle.  Experimental 
values  are  shown  on  Figure  5  for  an  actual  compressor  inlet .  At 
the  impeller  exit,  c  is  defined  as  shown  in  Appendix  A-2B.  A 
good  exit  profile  yields  97  per  cent,  a  poor  one  92  per  cent.  Between 
the  inlet  and  exit,  linear  or  near-linear  variation  is  assumed. 
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For  all  high  performance  and  high  Mach  number  compressors, 
the  blades  are  designed  as  thin  as  possible.  Manufacturing  con¬ 
siderations  establish  the  minimum  blade  thickness  at  the  shroud 
while  stress  and  vibration  considerations  establish  the  required 
blade  taper. 

The  density  in  the  above  equation  must  be  converged  upon  by 
iteration.  As  a  first  trail,  one-dimensional  lsentropic  densities 
can  be  used  for  the  inlet  portion  upstream  of  the  impeller;  inside 
the  impeller,  equation  A-8  is  used  in  conjunction  with  an  assumed 
variation  of  w  (in  general,  the  ultimate  desired  variation  is 
satisfactory).  L.tici  , as  the  design  oecoxnes  more  firmly 
established,  densities  and  velocities  will  be  recomputed  until 
satisfactory  accuracy  is  obtained. 

The  flow  picture  thus  found  is  used  as  a  preliminary  criterion 
of  design.  With  some  experience  in  the  interpretation  of  velocity 
distribution  some  modifications  can  be  anticipated,  discrimination 
can  be  made  between  basically  good  and  bad  designs  and,  finally, 
a  preliminary  selection  can  be  made  of  a  promising  design. 

As  indicated  earlier,  the  meridional  flow  field  obtained  thus 
far  is  an  approximate  one.  Corrections  must  be  made  for  the 
three-dimensional  nature  of  the  flow  and  for  the  effect  of  nonuniform 
work  addition.  Several  approximate  corrections  for  the  former 
effect  are  available;  one  such  correction  is  given  in  ref.  (1).  This 
conection  has  the  disadvantage  of  being  able  to  correct  the  wall 
velocities  only.  A  more  preferred  correction  is  given  ia  Appendix  C. 
Here  it  is  demonstrated  that  the  three-dimensional  incompressible 
flow  pattern  may  be  deduced  from  the  two-dimensional  incompressible 
pattern  by  taking  the  velocity  gradient  along  any  potent* al  line  of 
the  three-dimensional  flow  to  be  ?/4  o.'  the  corresponding  gradient 
of  the  two-dimensional  flow.  In  prat  tice,  this  correction  is  not 
large  and  its  effect  on  velocity  gradients  along  a  streamline  is 
negligible.  Care  must  be  taken,  ho*. ever .  to  apply  this  correction 
in  the  region  immediately  upstream  ci  the  impeller  leading  edge. 

The  correction  will  always  l>e  such  as  to  increase  the  flow  incidence 
near  the  shroud  by  an  amount  which  cannot  be  neglected. 
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A  sec  aid  correction  must  be  made  for  the  eL  <^'t  cf  non- 
uniform  work  addition  along  the  blade  span.  This  correction  is 
made  simultaneously  with  re-satisfaction  of  radial  equilibrium 
in  the  meridional  plane.  This  correction  is  necessary  because 
of  tbs  streamline  distortion  resulting  (1)  from  the  introduction  of 
the  three -dimensional  correction  and  (2)  from  the  effect  of  variable 
blade  blockage  from  hub  tc  shroud. 

The  equations  of  Appendix  B  express  the  condition  of  radial 
equ  Uibrium  for  the  case  of  transfer  of  mechanical  enei-gy  to  or 
from  an  incompressible  fluid.  Their  solution  requires  the 
knowledge  of  the  shape  o!  the  blade  mean  camber  line.  This 
solution  will  be  discussed  in  part  (b)  below.  The  solution  of  these 
equations  also  requires  the  knowledge  of  the  slope  and  radii  of 
curvature  of  the  meridional  streamlines.  This  Information  must 
be  obtained  by  an  iterative  process  in  which  successive  streamline 
plots  are  made  and  streamline  slopes  and  curvatures  obtained 
for  substitution  In  the  equations.  From  experience,  it  has  been 
found  that  five  or  more  iterations  are  necessary  for  convergence . 

(b)  Blade- to -Biaae  Solution 

All  calculations  up  to  this  point  have  been  carried  out  under 
the  assumption  that  the  flow  is  symnu  ♦rical  about  the  axis  erf 
rotation  i.e.,  that  flow  corresponding  t'  an  infinite  number  of  blades. 
Reference  19  indicates  that  this  flow  is  representative  of  the  mean 
flow  between  blades  for  an  impeller  with  a  finite  number  of  blades . 

A  separate  procedure  is  required  to  calculate  the  blade  loading  of 
an  impeller  with  a  finite  number  of  blades.  This  blade  ioading  is 
used  as  the  basic  criterion  of  impeller  design.  Appendix  D  gives 
equations  both  for  the  direct  and  the  inverse  problem.  The  direct 
solution  expresses  the  blade  loading  of  a  given  impeller  operating  with 
given  inlet  and  exit  conditions.  The  inverse  solution  enables  one 
to  calculate  the  blade  shape  required  to  produce  a  given  desired  loading. 
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The  blade -tc-blade  solution  must  oe  carried  out  prior  to  the 
corrections  for  vorticity  and  it  usually  most  be  repeated  afterward. 
Fortunately,  small  changes  m  the  blade  shape  have  but  little 
effects  upon  the  vorticity  correction,  so  that  repetition  of  this 
lengthy  correction  calculation  is  seldom  warranted.  Final  results 
are  usually  presented  as  relative  velocities  (mean  pressure,  and 
suction  side!  along  the  hut.,  shroud,  and  50  per  cent  streamlines. 

This  procedure  for  the  blade- to- olade  solution  is  based  on 
channel  theory  and  is,  therefore,  unreliable  in  the  regions  near 
the  leading  and  the  trailing  edges.  Near  the  leading  edges,  a 
separate  more  detailed  procedure  is  involved.  Near  the  trading 
edges,  corrections  have  to  be  introduced  to  account  for  the 
impeller  slip.  The  flow  angles  obtained  differ  from  the  corresponding 
blade  angles  principally  because  of  the  rotational  nature  of  the 
relative  flow  and  because  of  the  boundary  layer  thickness  on  the 
blade  suction  side.  Figure  6  is  a  summary  of  the  theoretical 
word  dene  at  NACA  for  mviscid  flow  between  rotating  radial  blades 
(ref.  171.  It  represents  the  theoretical  deviation  of  the  average 
flow  from  the  blade.  This  information  coupled  with  an  addrtional 
boundary  layer  deviation  can  be  used  *o  find  the  blade  shape  required 
to  produce  a  given  flow  angle  distribution.  The  total  deviation  at 
the  trailing  edge  is  known  (from  slip  factor  considerations)  and  using 
the  appropriate  curve  cf  Figure  6. the  deviation  due  to  the  bounuary 

layer  can  be  obtained  by  difference.  The  deviation  due  to  boundary 
iayer  is  assumed  to  decrease  linearly  irom  exit  to  inlet. 

Ttie  inducer  portion  very  near  the  leading  edges  must  ce 
designed  separately  to  establish  (1)  the  shape  of  the  leading  edge, 

(2)  the  blade  angle  variation  immediately  downstream  just  insiue 
the  impeller.  The  procedure  is  one  of  successive  approximations  in 
which  blade  cylindrical  developments  at  various  radii  must  be 
matched  together  aerodynamic  ally  and  structurally.  The  aerodynamic 
requirements  prescribe  the  variation  in  liK’idence  along  the  blade 
( »  2  degrees  in  the  supersonic  zone  and  .  7  degrees  ;n  the  low 
subsonic  zone*  according  to  the  data  of  references  25  to  33.  The 
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requirements  also  limit  the  amount  of  initial  deceleration  as 
the  flow  enters  the  cascade.  In  subsonic  cascades,  a  ratio  of 
cascade-throat  area  to  upstream  relative  flow  area  of  aboi  t 
1 .08  is  usually  found  best.  In  transonic  cascades  such  a  ratio 
is  not  definitely  established,  but  is  believed  to  be  somewhat  lower . 

The  structural  requirements  call  for  radial  blade  elements  and 
a  prescribed  blade  taper.  In  general,  all  these  conditions  tend 
tc  overdefine  the  design  of  the  blade  leading  edges,  in  which  case 
the  best  compromise  must  be  sought . 

A  new  approach  to  the  design  of  the  inducer  in  the  regions  of 
supersonic  flow  is  made  possible  from  the  present  study  of 
transonic  inducers.  The  results  are  incomplete  in  that  only  the 
two-dimensional  problem  has  been  wo  iked  out  numerically. 
Nevertheless,  considering  this  restriction,  the  study  provides 
means  for  finding  the  optimum  operating  incidence  for  a  given 
supersonic  inducer.  Conversely  the  blade  shape  required  to 
best  meet  fixed  upstream  conditions  can  be  calculated.  Also, 
some  indication  of  the  possibility  of  boundary  layer  separation 
as  a  result  of  shock  interaction  could  conceivably  be  deduced, 
although  the  task  would  be  laborious  and  the  results  of  doubtful 
validity . 

(c)  Inlet  Design 

The  description  of  the  flow  in  the  inlet  is  obtained  in  a  manner 
identical  to  that  described  earlier  for  the  impeller .  Starting  from 
the  inviscid,  ir  rotational ,  quasi-three  dimensional  ^quasi -compressible 
flow,  corrections  are  introduced  (1)  for  the  boundar^&splacement 
thickness  by  a  Pov  factor  varying  linearly  from  1.0  outside  to  0.93 
at  the  start  of  the  impeller,  and  (2)  for  the  three-dimensionality 
of  the  flow  as  previously  described.  With  high  meridional  Mach 
numbers,  the  transverse  gradients  in  velocity  are  corrected  for 
compressibility  by  using  the  Karman  Tsien  relation  referred  to  the 
one-dimensional  Mach  number .  This  correction  is  made  with  some 
reservation,  however,  insofar  as  there  is  some  doubt  of  its  validity 
for  certain  applications. 
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Finally,  graphical  matching  between  the  streamlines  in  the 
mlet  and  those  m  the  impeller  most  be  made.  A  good  Inlet  is 
considered  to  be  one  in  which  the  How  is  accelerating  everywhere 
and  in  which  leading  edge  velocity  gradients  are  held  down  as  much 
as  is  compatible  with  the  over-all  geometrical  requirements. 

2.  Flow  Analysis  by  Numerical  Procedures 

This  problem  consists  of  making  an  analysis  ol  the  flow  with 
given  geometry  by  direct  solution  of  the  basic  equations  of  fluid  flow. 
Reference  52  combines  the  basic  equations  of  fluid  How  into  a  convenient 
form  as  given  as  C-l  m  Appendix  C.  For  the  problem  at  hand,  equation 
C-l  may  be  simplified  to  equation  C-7  which  is  valid  for  compressible, 
lsentropic  flow  through  an  impeller  with  radial  elements.  Equation  C-l 
may  be  simplified  further  to  equation  C-3  in  the  inlet  portions. 

The  solubon  presently  under  way  attempts  to  solve  equations  C-3 
and  C-7  by  relaxation  techniques  on  high  speed  computers.  The  basic 
equations  are  expressed  in  finite  differences  and  solved  by  being  applied 
to  a  grid  of  abo«»»  600  points  in  the  impeller  passage.  In  order  to 
permit  computation  wuthin  the  capacity  of  the  IBM  650  machine,  a 
procedure  of  dual  iteration  is  used.  For  each  sequence  of  iteration 
(intended  to  satisfy  the  stream  function)  the  density  is  maintained  at 
the  value  calculated  from  the  previously  computed  velocities.  With 
new  values  of  the  stream  function,  new  densities  are  obtained,  and  so  on. 

Unusual  difficulty  is  introduced  at  the  impeller  leading  and  trailing 
edges  because  discontinuities  exist  there  both  in  the  direction  of  the  flow 
and  m  the  form  the  governing  equations.  This  condition  necessitates 
division  of  the  impeller  as  picturea  on  Figure  7.  each  section  to  be 
studied  separately.  Equation  C-3  for  compressible,  irrotational  flow  is 
solved  in  region  AA'B'B  and  the  results  are  considered  valid  in  the  region 
AA'C'C.  Equation  C-7  for  compressible  rotational  polytropic  flow  in 
the  blade  passages  is  solved  in  the  region  CC’T’T.  T.'e  ooundary  conditions 
along  CC’  are  known  from  the  above  calculations .  The  variation  of  flow 
direction  between  CC'  and  LL'  is  arbitrarily  prescribed  a  priori. 

Likewise  the  boundary  conditions  at  TT’  are  arbitrarily  prescribed.  The 
resuits  are  considered  valid  in  the  region  BB*D*D. 
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Region  CC’B'B  is  then  investigated.  The  boundary  conditions 
are  known  but  tne  problem  le  one  ts  which  the  condition  along  LL'  must 
be  studied  by  successive  approximations  until  both  sides  of  LL'  6how 
continuity  in  flow  velocity  (including  direction) . 

A  like  situation  exists  at  the  trailing  edge  where  the  two  regions 
EE'T’T  and  TT'D'D  must  be  matched. 

As  described,  it  is  obvious  that  this  complete  solution  is  an 
undertaking  of  considerable  magnitude,  and  is  too  time  consuming  for 
design  purposes . 

C .  SUMMARY  OF  DESIGN  PROCEDURE 

For  clarity,  the  over -all  design  procedure  may  be  summarised  as 
follows  in  the  order  in  which  it  must  be  carried  out . 

1 .  Impeller  External  Geometry 

(a)  Establishment  of  exit  velocity  triangle.  This  is  based  on: 

(1)  Selected  exit  blade  angle,  as  established  by  desired 
compressor  characteristics . 

(2)  Selected  impeller  degree  of  reaction. 

(3)  Estimated  slip  factor. 

(b)  Selection  of  mean  exit  diameter .  This  depends  upon  the 
compressor  application,  the  matching  of  the  turbine  and  compressor, 
the  over -all  limitation  in  size,  weight  etc. 

(c)  Proportioning  of  the  impeller  is  made  by  using  specific 
speed  and  minimum  inlet  relative  Mach  number  considerations  to 
find  the  hub/tip  ratio. 
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(d)  Minimum  inducer  bub  diameter  follow*  from  cutter 
clearance  requirements  for  the  estimated  number  of  blades. 

(e)  With  the  foregoing  established  the  exit  width  is  that 
required  to  pass  the  mass  flow . 

2.  Impeller  Internal  Geometry 

(a)  Field  plot  studies  are  made  of  a  series  of  potentially 
promising  meridional  shapes.  For  each  shape,  the  velocities 
are  calculated  from  the  streamline  or  potential  line  pattern 
and  are  subsequently  corrected  for 

(1)  Compressibility:  outside  the  impeller,  one- 
dimens  ionni  compressibility  corrections;  inside  the  impeller, 
densities  are  calculated  from  a  desired  variation  in 
relative  velocities. 

(2)  Blade  and  viscous  clogging  effects. 

(3)  Three-dimensionality  . 

(4)  Effects  of  nonuniform  work  addition .  At  this 
stage,  this  correction  must  be  estimated  and  its  value 
recalculated  later . 

(b>  The  shape  of  the  blade  is  calculated  by  a  direct  method 
to  provide  the  desired  loading  at  the  shroud. 

(c)  In  the  process  ol  the  design,  undesirable  features  of 
certain  meridional  shapes  may  become  apparent.  These  shapes 
are  eliminated.  More  refined  analyses  are  maoe  of  the  remaining 
shapes.  They  include:  two-dimensional  density  corrections  in 
the  inlet,  final  calculated  impeller  densities,  and  vorticity  effects. 

(df  Careful  design  of  leading  edges. 

(e)  The  entire  design  should  then  be  exactly  analysed  by  the 
numerical  method. 
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m  DESIGN  OF  IMPELLERS 

As  concluded  in  the  preceding  section,  Impeller  Design  Procedure, 
a  sufficient  number  of  Impellers  should  be  designed  to  investigate  the 
effects  of  three  primary  variables;  namely,  the  impeller  length,  the 
blade  loading  distribution,  and  the  exit  rone  angle.  A  total  of  ten 
impellers  have  been  designed.  They  are  subdivided  Into  three  classes, 

I,  n,  and  IV,  including  three  or  four  impellers  each.  Classes  I  and  n 
feature  impellers  with  60-degree  exit  cone  angle,  the  former  involving 
comparatively  long  Impellers  with  lighter  blade  loading,  the  latter  shorter 
impellers  with  greater  blade  loading.  Class  IV  features  impellers  wiih 
30-dcgree  exit  cone  angle  and  lighter  loadings.  Each  class  Includes  at 
ieaet  one  impeller  with  the  higher  loading  at  the  inducer  and  designated 
by  the  letter  "A";  one  impeller  nominally  with  uniform  loading,  without 
letter  designation:  and  one  impeller  with  the  higher  loading  displaced 
toward  the  discharge  and  designated  by  the  letter  "B*' .  The  main  over-ail 
features  of  each  impeller  are  tabulated  in  Table  I. 

A  summary  of  the  design  conditions,  assumptions,  ana  objectives  is 
given  in  Table  n. 

A.  Basic  Design  Features  Common  to  All  Configurations 

Inasmuch  as  the  present  program  is  concerned  with  the  optimization  of 
the  impeller  internal  geometry,  it  is  essential  that  all  impellers  be  sub¬ 
jected  to  identical  inlet  and  exit  conditions. 

(1)  Velocity  Triangles 

The  mean  or  one-dimensional  inlet  and  exit  velocity  triangles  are  ail 
identical.  Because  of  the  nature  of  the  loading  distribution,  however,  small 
variations  are  present  at  shroud  and  hub. 
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The  mean  inlet  and  exit  triangles  are  show;*  on  Figure  8 .  In 
accordance  with  the  criteria  formulated  under  Section  n,  the  exit 
blade  angle  lb  very  nearly  radial.  The  exit  Mach  number  is  slightly 
supersonic,  but  may  easlty  oe  reduced  to  a  safe  vilue  for  eventual 
matching  with  a  vaned  diffuser.  The  slip  factor  was  estimated  to  be 
0.81. 


The  mean  Inlet  diagram  corresponds  to  the  local  velocity  triangle  at 
the  50  per  cent  streamline.  This  designation  Indicates  that  over  50  per  cent 
of  the  inlet  flow  Is  supersonic.  A  typical  inlet  relative  Mach  number 
variation  with  radius  is  as  shown  In  Figure  9,  while  Fi  gure  10  shows  a 
similar  variation  in  inlet  absolute  Mach  number. 

A  relative  mean  deceleration  ratio  Wl  =  0.625  wai  selected.  This 

w  i 

represents  a  considerable  amount  of  deceleration,  especially  in  view 
of  the  fact  that  the  deceleration  ratio  at  the  shroud  will  be  considerably 
lower  and  still  lower  along  the  blade  suction  surface.  In  view  of  the 
design  philosophy  outlined  in  Section  11,  however,  a  value  of  reaction 
slightly  above  50  per  cent  is  desirable,  and  the  present  reaction  Is 
59  per  cent.  It  is  felt,  furthermore,  that,  if  at  a  later  date  it  becomes 
apparent  that  excessive  deceleration  is  present,  this  variable  may  be 
investigated  by  successive  decreases  of  the  exit  width.  This  procedure 
thus  considers  eventual  matching  of  this  impeller  with  a  diffuser. 

(2)  Lxit  Cone  Angle 

For  impeller  Classes  I  and  Q,  a  60-degree  exit  cone  angle  was 
selected.  This  angle  is  felt  to  be  best  for  the  present  specific  speed.  The 
effect  of  smaller  angles  upon  performance  cannot  be  fully  evaluated 
beforehand  but  will  be  investigated  through  impeller  Series  IV.  A  smaller 
angle  will,  however,  obviously  result  in  a  heavier  wheel. 
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(3)  Diameter  Ratios 

As  dictated  by  the  consideration  cf  specific  speed,  the  ratio  of 
Inducer  tip  diameter  to  mean  impeller  tip  diameter  is  0.713  for  all 
impellers.  The  inducer  hub/tip  diameter  ratio  was  designed  as  low  as 
compatible  with  blade  root  thickness  and  cutter  dimensions:  it  is  0.322. 

The  impeller  exit  hub/tlp  diameter  ratio  is  set  by  the  desired 
distribution  of  work  between  the  hub  and  shroud.  Because  otf  the  loss 
distribution  inside  the  impeller  (the  majority  of  which  is  estimated  to  be 
at  the  shroud  because  of  supersonic  relative  speeds),  high  ratios  of 
deceleration,  blade  leakage  effect,  and  accumulation  of  secondary  flow, 
more  work  must  be  done  by  the  impeller  along  the  blade  tip,  if  a  uniform 
distribution  of  total  pressure  is  to  be  achieved  at  the  exit.  This  work 
increase  at  the  tip  was  set  at  three  per  cent  more  than  the  mean  and  at 
the  root  three  per  cent  less  than  the  mean  which  brought  about  a  hub/tip 
ratio  at  the  exit  of  0.94,  this  ratio  being  the  same  for  all  impellers. 

(4)  Inlet  Arrangement 

Hie  shape  of  the  inlet  duct  to  the  impeller  has  been  selected  as  one 
turning  the  flow  from  a  radial  direction.  Some  controversy  exists  as  to 
the  relative  merits  of  designing  for  axial  or  radial  inflow.  Clearly,  an 
axial  inlet  is  superior  from  the  standpoint  of  standardization  and 
mechanical  simplicity.  On  the  other  hand,  the  radial  inlet  is  more 
representative  of  the  majority  of  small  gas  turbine  installations .  Further¬ 
more,  mixed-flow  compressors  require  a  straddle- bearing  arrangement; 
an  axial  inlet  would  lose  part  of  its  advantage  of  simplicity  in  that  it  would 
have  to  be  a  structural  element  and  provide  front  bearing  space  and  support. 

It  thus  appears  justified  to  use  a  radial  inlet  for  the  present  research  program. 
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B.  Internal  Aerodynamic  Design  of  Impellers 

In  accordance  with  the  conclusions  of  Section  I  ,  the  problem  at  hand 
is  to  produce  for  each  class  of  impeller  a  series  of  designs  which  will 
permit  variation  of  the  distribution  of  blade  loading  by  varying  only  one 
physical  variable  (meridional  shape  or  Llade  camber)  and,  yet,  will  be 
near  optimum  throughout.  From  a  study  made  to  investigate  which  of 
the  two  to  vary  it  was  found  most  desirable  to  hold  the  meridional  shape 
constant  and  change  the  blade  camber  distribution.  For  each  impeller 
class,  four  to  five  meridional  shapes  were  studied  until  one  was  found 
optimum  with  a  blade  shape  producing  uniform  loading  while  near  optimum 
for  otner  loading  distributions.  This  shape  was  retained  for  the  final 
design  and  subjected  to  the  more  refined  design  procedure  for  each 
particular  desired  blade  loading.  The  details  of  the  procedure  are  fully 
described  in  Section  n. 

Some  additional  work  was  expended  upon  the  design  ot  the  impeller 
leading  edges.  Since  transonic  impellers  have  the  undesirable  characteristic 
of  being  capable  cf  operation  only  in  a  very  narrow  flow  range,  a  more  refined 
design  approach  must  be  made  at  the  inlet,  if  design  flow  is  to  be  met. 

This  is  particularly  so,  since  there  has  been  evidence  that  the  present 
design  procedure  tends  to  overestimate  the  inlet  shroud  velocities  somewhat. 
Among  other  things,  a  study  was  made  to  investigate  me  effect  oi  successive 
readjustments  of  the  streamline  radii  in  the  meridional  plane,  since  this  is 
not  always  done  in  the  basic  design  procedure.  The  resultant  shroud 
velocities  were  found  to  be  reduced  by  as  much  as  10  per  cent,  compared  to 
the  originally  calculated  values . 

From  these  results,  it  was  felt  necessary  to  undertake  a  more  exact 
solution  of  inlet  flow.  This  effort  involves  the  numerical  solution  of  the 
equations  of  motion  on  an  IBM  650  high  speed  computer .  This  work  is 
being  done  in  two  phases,  in  increasing  degrees  of  difficulty.  The  first 
phase  treats  of  the  three-dimensional,  ir rotational  flow  in  the  inlet  and 
impeller  passage  without  blades  •  itti  o  VTGI  L  is  complete.  The  second  phase 
treats  of  the  same  problem  but  for  compressible  flow.  Contrary  to  what  was 
originally  expected,  the  second  of  these  phases  has  not  yet  progressed  far 
enough  for  results  to  be  incorporated  herein.  While  these  results  are  not  as 
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yet  complete,  they  afford  some  information  in  the  form  of  predicted 
velocity  upstream  of  the  leading  edge .  The  combined  effect  of  the  iteration 
of  the  streamlines  and  the  more  exact  solution  of  the  three  dimensional 
flow  indicates  a  shroud  meridional  velocity  about  10  per  cent  lower  than 
originally  calculated.  The  net  result  is  that  a  forward  extension  of  the 
leading  edge  is  required  for  the  design  flow  to  meet  the  blades  at  the 
correct  incidence.  The  design  incidence  was  allowed  to  vary  from  3 
degrees  at  the  shroud  to  5  degrees  at  the  relative  sonic  point  to  7  degrees 
at  the  hub.  This  incidence  variation  is  the  same  for  all  impellers.  This 
effect,  together  with  the  variation  of  inducer  blade  loading  among  impellers, 
causes  small  differences  in  the  amounts  ad  inducer  extension  forward  of 
the  original  position.  Figure  11  illustrates  typical  shapes  of  leading  edges. 

It  must  be  emphasized  that  these  modifications  will  not  detract  from 
the  objectives  of  the  original  designs.  True,  the  reported  velocities  may 
differ  by  as  much  as  10  per  cent  from  those  calculated  by  the  more  exac' 
solution,  but  the  gradients  and,  above  ail,  the  loading  of  each  impeller 
relative  to  ary  other  will  be  unchanged. 

More  will  be  said  about  this  flew  analysis  in  the  final  report  when  ail 
results  are  available.  In  the  meanwhile,  Appendix  C  summarizes  the 
method  of  attack  and  produces  the  pertinent  equations. 

C .  Geometry  of  Impellers 

Description  of  the  geometry  of  the  impellers  involves  blade  coordinates 
and  meridional  shape  coordinates.  Both  are  presented  in  accurate  tahiilar 
and  graphical  forms  against  the  axial  coordinate  Z .  A  "smoothing  out"  of 
the  tabulated  coordinates  is  necessary  before  machining  of  impellers  can  be 
started.  That  is,  the  coordinates  must  be  refined  to  a  point  where  the  first 
and  second  derivatives  are  smooth  and  continuous  functions.  Only  the 
coordinates  for  Impellers  I,  IB,  and  II,  which  are  to  be  fabricated  first  have 
been  so  smoothed  out.  The  meridional  coordinates  ofthe  inlet  and  Impeller 
Series  I  and  11  are  tabulated  on  Tables  IQ  through  V  Drawings  of  the 
meridional  shapes  are  seen  as  Figures  12  to  14.  Tabulation  of  "smoothed 
out”  blade  coordinates  are  given  in  Tables  VI  through  VQI  and  plots  for  the 
blade  angular  coordinates  of  all  impellers  are  shown  on  Figures  15  through  24. 
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#D.  Impeller  Aerodynamic  Characteristics 

Figures  25  through  64  picture  the  impeller  cnaracteristics  fcr  all  designs. 
They  are  presented  in  the  form  of  blade  loading  curves  expressed  as  relative 
velocities  and  In  terms  of  mean  pressure  along  the  shroud  and  along  the  hub . 
All  curves  are  typical'of  impellers  with  radial  blade  elements.  Since  the 
entire  blade  shape  is  fixed  when  the  blade  contour  is  prescribed,  say  at  the 
shroud,  there  is  not  the  freedom  to  control  the  loading  both  at  the  shroud 
and  hub.  The  desired  loading  is  built  in  at  the  shroud,  since  this  region  is 
the  most  influential .  Appreciable  hub  loading  in  the  inducer  cannot  be 
achieved  without  overloading  the  shroud. 

(lj  Impeller  Series  1 

Of  all  three  meridional  configurations  that  of  Impeller  Series  I  proved, 
during  design,  to  be  the  most  tracLible .  It  was  found  that  the  desired  loadings 
could  be  achieved  comparatively  easily  .  li  is  felt  that  this  design  is  the  most 
conservative  and  the  one  best  suited  to  the  design  specific  speed. 

(2)  Impeller  Series  n 

With  this  series,  the  impeller  design  is  short  and  crowded  together 
resulting  in  a  certain  lack  of  freedom  and  flexibility .  The  uispla/ed  loadings 
were  obtained  only  at  the  expense  of  designing  with  thinner  blades.  They 
are  so  thin,  in  fact,  that  serious  machining  problems  are  anticipated.  All 
data  given  herein  are  for  the  design  with  thin  blades,  but  the  possibility 
remains  that  a  slight  increase  in  thickness  may  have  to  be  allowed,  if 
difficulties  are  encountered  in  fabrication.  This  impeller  was  designed  with 
18  blades  nstead  of  17  for  Impeller  Series  I  because  of  the  somewhat  higher 
blade  loadings. 
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(3)  Impeller  Series  IV 

Changing  from  60  to  30  degrees  exit  cone  angle  was  originally  expected 
to  bring  about  a  shift  in  the  loading  distribution  relationship  between  shroud 
and  hub.  This  shift  proved  to  be  email  and  essentially  equivalent  hub 
loadings  were  achieved.  The  30-degree  exit  angle,  however,  is  too  low 
for  this  specific  speed.  This  was  apparent  in  that  it  was  found  difficult 
to  arrive  at  a  good  meridional  shape.  An  axially  short  design  called  for 
a  certain  amount  of  canting  at  the  discharge  as  shown  in  Figure  65.  This 
canting  produced  velocity  gradients  In  the  exit  region,  tending  to  unload 
the  blade  tip  and  causing  turbining.  In  order  to  avoid  this  condition,  a  longer 
impeller  was  needed  which  is  undesirable  from  a  weight  standpoint.  However, 
with  a  longer  impeller  loadings  are  lower  permitting  one  less  blade  than 
the  Series  I  impeller . 

A  common  inlet  passage  matches  all  configurations  and  is  designed  to 
provide  accelerating  flow  in  all  regions. 

E .  Impeller  Structural  Considerations 
(1)  Centrifugal  Stresses 

Analysis  of  the  stress  level  in  the  blades  and  in  the  hub  were  made  for 
a  steel  wheel  operating  at  10  per  cent  overspeed.  The  calculated  stresses 
are  those  due  to  rotation  alone  and  include  neither  vibratory  nor  thermal 
stresses.  Ttiese  latter  stresses  must  be  kept  in  mind  even  though  their 
magnitude  cannot  readily  be  calculated.  The  calculated  blade  stress  for 
Impeller  I  is  shown  on  Figure  66.  The  maximum  is  above  75, 000  psi  and 
occurs  at  the  root  at  Z  =  0.6.  The  critical  hub  stress  condition  is  shown 
on  Figure  67.  At  the  center  of  the  hub  the  maximum  stress  is  84, 000  p6i. 

The  corresponding  data  for  Impeller  Seriee  II  is  Indicated  on  Figure  SS 
and  69.  Similarly  the  maximum  stresses  are,  respectively,  86,000  and 
87,000  psi. 


WADC  TR56-589 


30 


CONFIDENTIAL 


CONFIDENTIAL 


Other  impellers  of  the  I  ar.d  II  Series  show  appreciable  differences 
in  their  relative  chord  length  of  the  blades,  thus  ezhibuiug  differences  in 
hub  stresses.  Hie  maximum  stress  point,  however,  i. 3  in  the  aft  portion  of 
the  impeller  where  the  blades  are  of  nearly  equal  length.  It  may  then  be 
stated  with  confidence  that  no  higher  maximum  stresses  will  be  found  in 
the  other  impellers. 

No  stress  analysis  has  yet  been  made  for  Impeller  IV .  Most  likely  it 
wiD  be  of  the  same  order  of  magnitude  or  lower . 

(2)  Critical  Speeds 

Calculations  were  made  for  the  impeller  critical  speeds  for  the  con¬ 
figuration  of  one  rear  (at  impeller  exit)  journal  bearing  and  a  ball  bearing 
in  the  front.  In  view  of  the  great  rigidity  of  the  impel*  er  and  short  stub 
shafts,  the  calculated  critical  speed  depends  largely  upon  the  estimated 
rigidity  of  the  journal  bearings.  With  the  best  available  estimate  the  first 
critical  mode  was  found  to  Le  10,400  rpm  and  the  second  mode  96,000  rpm. 
If  the  rigidity  weie  ten  times  higher  than  the  original  estimate,  the  first 
mode  would  only  be  27,  500  rpm.  It  appears  then  that  operating  speed  is 
sufficiently  far  away  from  any  critical  speed. 

(3j  Impeller  Thrust 

At  standard  conditions,  calculations  indicate  an  impeller  thrust  of 
1033  lb.  This  will  require  the  use  of  tandem  bearings  at  the  front  of  the 
impeller. 
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IV  STUDY  OF  TRANSONIC  INDUCER 

One  of  the  most  convenient  characteristics  of  low  Mach  number  fluid 
flow  (both  incompressible  and  compressible)  rests  in  its  ability  to  conform 
to  lees  than  optimum  geometrical  configurations.  Such  flows  will  tolerate 
a  comparatively  wida  variation  in  Incidence  upstream  of  a  cascade  with  only 
slight  effects  on  efficiency .  This  is  not  the  case  in  a  transonic  Inducer. 
Angles  and  effective  flow  variations  must  be  met  exactly,  or  sizeable 
penalty  in  efficiency  will  be  incurred.  As  a  result,  conventional  design 
methods  are  no  longer  adequate;  more  exact  approaches  for  predicting 
flow  behavior  must  be  devised.  In  high  specific  speed  machines  such  as 
axial  compressors,  it  is  quite  evident  that  losses  in  the  inlet  will  have  a 
large  effect  on  the  07er-all  efficiency,  since  the  inlet  relative  dynamic 
head  is  a  large  fraction  of  the  total  head  of  the  compressor.  In  lower 
specific  speed  units,  tests  at  AiResearch  on  high  Mach  number  machines 
of  radial  and  mixed-flow  types  have  shown  that  unsatisfactory  distribution 
c!  incidence  along  the  leading  edge  is  apt  to  disrupt  completely  the  flow  inside 
the  impeller  as  evidenced  by  tne  distorted  shape  of  the  velocity  prefile  out 
of  the  impeller.  It  is,  therefore,  of  utmost  importance  that  work  be  done 
toward  improving  the  general  understanding  of  the  flow  phenomena  in  a 
transonic  cascade.  Only  then  will  it  become  possible  to  design  blading  such 
that  diffusion  takes  place  without  separation  or  such  that  shock  patterns 
can  be  accommodated  witnout  extraneous  losses. 

A.  Review  of  Previous  Work 

Up  to  this  time  (as  indicated  by  available  publications)  the  work  done 
toward  providing  this  information  may  be  subdivided  as  follows: 

1 .  Mathematical  formulation  of  the  basic  equations  describing 
the  fluid  motion  between  blades  of  an  impeller  with  a  finite  number 
of  blades  have  been  made.  Reference  43  deals  independently  with  two 
relative  stream  surfaces,  one  radial  and  the  other  t^uigentlal  upstream 
of  the  impeller .  Solution  (leading  to  a  description  of  the  flow)  requires 
a  lengthy  iterative  procedure  until  the  basic  equations  are  everywhere 
satisfied  for  both  surfaces.  No  complete  solution  is  known  to  have 
been  carried  out. 
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2.  For  the  two-dimensional  case,  several  qualitative  or 
quantitative  solutions  have  been  worked  out,  mainly  by  graphical  methods 
and  lor  supersonic  Mach  numbers  of  the  order  of  2G.  References  3% and 
44  are  typical. 

3 .  Considerable  experimental  work  has  been  done  (as  In  references 
44  and  45)  on  supersonic  cascades  by  observation  of  the  flow  with  optical 
means.  On  rotating  cascades,  successful  Sehlieren  pictures  have  been 
taken  by  NACA  fc»*  a  compressor  rotor  with  supersonic  flow  at  the  tip. 

Hie  problem  at  hand  differs  from  otner  similar  ones  treated  In  the 
literature  in  that  (1)  the  blade  solidity  is  high  compared  to  that  of  a  system 
of  axial  /anes,  (2)  the  relative  Mach  number  range  up  and  down  the  leading 
edge  is  large,  from  low  supersonic  at  shroud  to  low  subsonic  at  hub,  and 
(3)  the  superbonic  tip  Mach  number  is  lower,  closer  to  1.0,  in  general,  than 
in  other  cascade  «o  investigated. 


B.  General  Procedure 

The  following  general  procedure  here  prescribed  involves  a  back -and -forth 
iterative  procedure  on  two  families  of  surfaces  one  of  them  being  a  series  of 
surfaces  of  revolution,  on  which  the  blade -to -blade  flow  must  be  satisfied, 
the  other  being  a  series  of  surfaces  generated  by  a  radial  element  in  the 
fluid  moving  parallel  to  the  blade  camber  surface  on  which  radial  equilibrium 
must  be  satisfied.  The  blade -to -blade  flow  is  to  be  solved  at  four  different 
radii,  two  radii  where  the  relative  flow  is  supersonic  and  two  where  it  is 
subsonic .  In  the  supersonic  region,  graphical  solution  will  be  used  by  an 
adaptation  of  the  conventional  field  buiution  while  in  the  subsonic  region  an  analog 
field  plotter  will  be  used  and  the  results  subsequently  corrected  for  com¬ 
pressibility  and  three-dimensional  effects.  Radial  equilibrium  is  to  be  satisfied  at 
three  angular  positions,  blade  suction,  and  pressure  sides,  and  mid-channel. 

The  procedure  is  similar  to  that  used  in  the  solution  of  the  flow  equations  for 
the  axi-symmetiic  condition.  Its  objective  is  to  find  the  radial  variation  in 
the  thickness  of  the  stream  sheets  adjacent  to  the  four  selected  surfaces  of 
revolution.  This  thickness  is  expected  to  vary  in  he  tangential  as  well  as  in 
the  axial-radial  direction. 
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The  Iterative  procedure  consists  a i  recalculating  the  blade-to-blade 
flow  allowing  for  the  thickness  variation  in  the  third  dimension  and  then 
repeating  the  radial  flow  equilibrium  calculation  to  obtain  new  stream 
sheet:  thickness  variation.  Repetition  c*  this  procedure  must  be  continued 
until  convergence  is  met  everywhere . 

Boundary  conditions  for  upstream  and  downstream  are  established  by 
an  independent  design  analysis.  Upstream  in  the  inlet  axi-symmstrlcal 
flow  applies  while  downstream  the  blade -to -blade  flow  is  computed  by 
channel  theory .  The  downstream  boundary  must  be  taken  only  where  the 
relative  flow  is  subsonic,  since  the  channel  flow  analysis  does  not  distinguish 
between  subsonic  and  supersonic  flows . 

C .  Extent  of  Work  Completed 

The  original  general  approach  to  this  part  of  the  over-all  program 
was  to  carry  on  this  inducer  work  simultaneously  with  the  Impeller  design 
study  but  on  a  lower  priority  basis .  As  work  progressed,  it  soon  became 
apparent  that  the  inducer  work  was  a  problem  requiring  a  formidable  effort 
to  bring  it  to  completion.  It  also  became  obvious  that  only  an  enormously 
simplified  version  of  the  above  procedure  would  have  usefulness,  if  it  were 
to  be  incorporated  as  part  of  the  design  procedure .  It  was  therefore  decided 
to  carry  this  work  only  as  far  as  time  would  allow  and,  rather  than  attempt 
to  reach  a  complete  detailed  numerical  solution,  look  for  possible  features 
useful  for  design  and  tor  a  better  general  understanding  of  the  phenomenon. 

The  Initial  assumption  for  the  first  attempt  at  the  blade-to-hlade 
solution  was  taken  as  a  constant  stream  sheet  thickness .  Even  this  two- 
dimensional  solution  requires  a  compound  iterative  procedure  in  the 
supersonic  zone.  It  was  carried  out  for  cne  stream  surface  only,  the  one  at 
the  shroud.  As  anticipated,  with  the  imposed  condition  of  two -dimensionality, 
it  was  not  possible  to  satisfy  the  downstream  conditions  completely ,  that  is, 
the  three-dimensional  effects  must  be  sizeable.  In  spite  of  this  result,  some 
understanding  was  gained  as  to  the  general  nature  of  the  flow  pattern .  Means 
of  computing  the  optimum  blade  incidences  have  been  found,  and  some  insight 
has  been  gained  about  the  relationships  existing  between  blade  shapes  and 
the  form  of  the  shock  patterns. 
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D.  Significant  Qualitative  Results 

Two  possible  basic  shock  configurations  were  found  to  exist;  one  is 
represented  by  a  line  such  as  LCR  on  Figure  70  (curved-oblique  shock 
configuration!,  the  other  by  a  line  such  as  LCR' (curved-normal  shock 
configuration).  Both  the  shape  of  the  shock  (whether  normal  or  oblique) 
and  its  location  upstream  of  the  leading  edge  can  be  calculated.  The  shape 
and  location  of  the  shock  are  functions  of  the  blade  shape,  the  upstream  flow 
oSach  number  and  direction,  and  the  downstream  conditions . 

1 .  The  shape  and  location  of  the  left-hand  branch  of  the 
shock,  for  a  given  cascade  under  given  upstream  conditions,  is 
entirely  determined  by  the  requirement  of  periodicity  upstream 
of  the  blades.  The  calculation  can  be  made  independent  of  the 
downstream  conditions  and  independent  of  the  shape  of  the  right-hand 
branch . 

2 .  The  shape  of  the  right-hai^d  branch  is  a  function  of  the 
downstream  conditions  and  requires  the  solution  of  the  entire  flow 
field. 

3.  Optimum  incidence  or  optimum  blade  shape  may  be  defined 
as  that  whx«.h  will  produce  at  R  or  R1  a  pressure  rise  least  likely 

to  bring  about  boundary  layer  separation.  Calculations  show  the 
reflected  curved-oblique  configuration  to  be  more  desirable  than  the 
normal  shock  configuration .  For  a  rigorous  application  of  this  criterion, 
a  complete  solution  of  the  flow  is  required.  With  lower  downstream 
pressures,  the  curved -oblique  shock  will  exist.  As  the  back  pressure 
is  raised,  the  shock  will  jump  into  the  normal  position.  When  this 
takes  place,  breakdown  of  the  flow  and  collapse  of  impeller  performance 
is  likely  to  take  place.  The  design,  then,  should  be  made  for  a  curved- 
oblique  shock  pattern . 
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4 .  To  the  two  possible  shock  configurations,  there  corresponds 
two  separate  mechanisms  by  which  the  flow  becomes  subsonic  for  an 
lnviscld  fluid.  When  the  shock  is  curved-normal,  such  as  LCR', 
everything  downstream  becomes  and  remains  subsonic .  When  the  shock 
is  curved-oblique,  such  as  LCR,  the  transition  from  super  to  subsonic 
takes  place  in  steps  as  sketched  on  Figure  70,  subsonic  flow  gradually 
filling  the  channel. 

It  is  possible  to  shape  the  blades  in  such  a  way  as  to  promote  the 
desirable  curved- oblique  configuration.  For  this  configuration,  the 
blade  curvature  must  be  such  that  the  oblique  shock  can  be  reflected. 
Insufficient  curvature  will  result  in  too  low  a  Llach  number  at  R  to  permit 
reflection  and  the  shock  will  jump  out  to  a  position  such  as  CR* . 
Excessive  curvature  will  slow  down  the  deceleration  process  and  result 
in  an  extensive  zone  of  supersonic  flow  between  the  blades . 

E.  Method  of  Analysis  of  the  Two-Dimensional  Supersonic  Flow 

The  procedure  of  analysis  outlined  below  is  applicable  to  the  two- 
dimensional  invisciil  flow  on  a  cylindrical  surface  of  large  radi  os.  The 
analysis  is  based  on  the  assumption  that  a  system  of  detached  shocks  is 
present  in  front  of  the  cascade.  Only  this  flow  pattern  is  logical  for 
flows  at  low  supersonic  llach  numbers .  Furthermore  such  a  pattern  has 
been  confirmed  experimentally  by  optical  observation.  The  shape  of  each 
detached  shock  differs  from  that  of  the  more  familiar  variety  existing  ahead 
of  a  single  airfoil  because  of  the  presence  of  adjacent  blades  and  because 
the  downstream  conditions  may  be  varied  independently  of  the  upstream 
conditions .  It  is  necessary  to  distinguish  between  the  center  portion,  the 
right  and  left  branches  of  the  shock.  Their  shapes  are  governed  by  separate 
considerations. 
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1 .  For  a  given  geometry  and  upstream  condition,  the  shape  of  the 
left  branch  LC  on  Figure  70,  and  its  position  ahead  of  the  leading  edge, 
is  dictated  purely  by  requirements  of  periodicity.  Thus,  at  poLit  P,  the 
Mach  number  must  go  through  a  discontinuity  from  the  value  it  has  along 
wave  W(  to  the  value  it  must  have  along  wave  W2  (which  must  be  the  same 
as  W8\)  If  the  local  Mach  numbers  were  known,  the  direction  of  the  shock 
at  F  (or  any  other  such  pulnt)  would  also  be  known . 

Assignment  of  the  numerical  values  of  Mach  numbers  a.t  various 
points  must  be  made  by  successive  approximations.  Assignment  of  a 
Mach  number  to  a  point  such  as  Q,  permits  the  drawing  of  all  expansion 
waves  W0,  Wt,  etc.  over  the  entire  field,  each  labelled  with  its  particular 
Mach  number.  At  all  points  such  as  P,  little  elements  of  shocks  can  be 
drawn  as  explained  above.  If  the  center  portion  of  the  shock  location  is 
selected,  the  entire  left  branch  can  be  sketched  in.  One  particular  Mach 
line  will  never  intersect  any  shock.  The  stream  velocity  and  direction 
along  this  line  must  be  equal  to  the  stream  Mach  number  and  direction 
far  upstream .  If  the  Mach  number  condition  is  not  met,  another  Mach 
number  value  18  assigned  at  point  Q  If  the  directional  requirement  is 
not  met,  the  position  of  the  shock  center  portion  is  moved  a  little 
upstream  or  downstream  relative  to  the  blade  leadaig  edge.  Thus, 
satisfaction  of  the  blade  geometry  and  the  upstream  conditions  fix  the 
shape  and  location  of  the  ieft  branch  of  the  shock. 

2.  The  shape  of  the  shock  right-hand  branch  (branch  CR)  must  be 
such  as  io  satisfy  the  downstream  conditions  previously  established  by 
channel  theory  analysis.  The  problem  is  one  of  calculating  a  flow 
satisfying  given  boundary  conditions. 

For  purposes  of  engineer Lig  computations,  it  is  perhaps  best  as  an 
initial  working  procedure  to  analyze  the  flow  by  graphical  means  rather 
than  attempt,  from  the  start,  to  set  up  a  complete  mathematical  solution. 
The  merits  of  this  ?«pproach  were  quit*  evident  in  the  calculation  of  the 
flow  on  both  sides  c.4  the  shock  ieft  branch  where  a  solution  could  be 
converged  upon  comparatively  rapidly.  The  flow  on  the  right  side  is 
mixed  subsonic  and  supersonic.  The  supersonic  regions  will  be  handled 
by  ihe  conventional  field  method  of  shocks  and  characteristic  lines  of 
finite  strength.  The  subsonic  regions  w;il  be  handled  bj  plots  of  stream¬ 
lines  such  that  continuity  and  radial  equilibrium  are  satisfied. 
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Solution  of  the  pattern  involves  a  method  of  successive  appro  delations 
with  adjustment  of  the  sonic  point  S  and  the  sonic  line  SN .  Adjustment 
of  SN  (or  BAIT)  is  required  to  satisfy  the  basic  flow  equations.  Changing 
the  location  of  S  must  be  made  until  the  desired  downstream  subsonic 
channel  conditions  are  met.  For  a  given  selection  of  point  8  and  line 
3N,  the  entire  region  SRN  cm  be  computed  in  the  conventional  way 
in  terms  of  the  properties  within  finite  patches,  between  shocks  and 
Mach  lines.  The  flow  picture  m  the  subsonic  region  is  necessarily  less 
precise  in  that  the  streamline  shape  must  be  calculated  and  plotted 
point  by  point.  The  reflected  shock  RN  must  terminate  as  a  normal 
shock  at  N,  s  nee  the  local  Mach  number  will  be  near  1.0.  Region  IfM’T 
is  computed  hi  a  similar  manner .  Whether  another  oblique  reflection 
takes  place  or  whether  the  supersonic  flow  finally  ends  in  a  normal  shock 
depends  on  the  downstream  condition.  At  the  point  where  the  /low  la 
subsonic,  the  velocities  on  both  the  suction  and  pressure  sides  cf  the 
blade  must  check  those  calculated  by  channel  theory.  If  this  is  not 
realized,  a  fresh  start  is  made  for  the  position  of  points . 

3 .  Little  concern  need  be  had  about  the  shock  center  region .  Its 
position  forward  of  the  leading  edge  follows  automatically  from  part 
(1)  above.  The  limits  of  the  subsonic  zone  and  the  shock  direction  at 
points  C  behind  the  shock  follows  from  parts  (1)  and  (2).  The  exact 
shape  of  that  portion  of  the  shock  has  little  effect  on  the  over-all  flow 
picture . 

F.  Calculated  Flow  Pattern 

Study  of  the  blade  loading  within  the  impeller  as  calculated  by  the  method 
of  Section  n  indicates  that  the  flow  remains  supersonic  on  the  blade  suction 
side  for  a  considerable  distance.  Thus,  if  downstream  conditions  are  to  be 
met,  a  system  of  curved-oblique  shocks  must  be  present.  Figure  71  pictures 
such  a  flow,  as  far  as  it  was  calculated,  and  is  the  result  of  a  series  of  iterations. 
The  flow  was  not  carried  all  the  way  downstream  when  it  became  apparent  that 
continuity  could  not  be  satisfied  in  the  subsonic  zone,  if  the  requirements  of 
two-dimensionally  alone  are  adhered  to.  This  indicates  that: 


WADC  TR56-589 


36 


CONFIDENTIAL 


CONFIDENTIAL 


1 .  A  two-dimensional  flow  is  not  possible  with  the  imposed 
upstream  and  downstream  conditions . 

2.  A  two-dimensional  solution  exists,  but  It  does  not  satisfy  the 
downstream  conditions  Such  a  solution  requires  a  frontal  normal 
shock  ruch  as  CR'  on  Figure  70.  This  solution  is  of  little  interest  here. 

3 .  An  illustrative  two-dimonslonal  flow  may  be  conceived  in 
vhich  its  upstream  and  downstream  conditions  are  met  approximately 
but  in  which  the  pressure  side  boundary  is  so  modified  as  to  permit 
satisfaction  of  continuity.  The  flow  pattern  of  Figure  71  conforms  to 
this  description.  The  loLl owing  interesting  points  were  encountered: 

(a)  The  shock  location  satisfying  beth  upstream  Mach  cumber 
and  direction  is  one  which  has  its  center  portion  very  near  the  blade 
leading  edge. 

(b)  Because  of  the  multiple  shocks  and  expansions  upstream  of 
the  leading  edges,  the  Mach  number  immediately  upstream  of  the 
center  portion  of  the  shock  Is  considerably  higher  than  In  the 
undisturbed  stream  (1.50  compared  to  1.32). 

(c)  The  possibility  of  the  right-hand  shock  branch  being  of  the 
curved  oblique  form  but  becoming  normal  at  the  section  surface  of 
the  adjacent  blade  was  Investigated  and  found  impossible. 

(d)  For  regular  shock  reflection  on  the  blade  suction  surface, 
the  local  incoming  Mach  number  must  be  above  1.5.  Lower  local 
Mach  numbers  will  cause  the  shock  pattern  to  alter  itself  to  the 
curved  normal  configuration. 

(e)  Regardless  of  the  position  of  the  sonic  point  s  and  the  shape 
of  the  streamline  sd,  continuity  cannot  be  satisfied  between  stream¬ 
lines  OO  *c  and  OO-d.  This  points  to  the  extreme  significance  of 
the  three -limcrsional  efiects. 
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G.  Process  of  Diffusion  in  Transonic  Impeller 

At  this  time,  it  is  of  great  intereet  to  compare  the  qualitlve  results 
of  the  foregoing  analysis  with  some  observed  behavior  of  AlResearch 
transonic  impellers.  Figure  72  shows  an  approximate  sketch  of  experimental 
shroud  pressure  data  for  a  typical  set  of  runs  of  a  transonic  radial  impeller. 

Of  particular  interest  is  the  fact  that  at  high  mass  flows,  extremely  low 
pressures  are  reached,  Indicating  very  high  local  Mach  numbers.  The 
humps  on  the  curve  indicate  the  presence  of  several  successive  shocks, 
separated  by  re-expansions.  As  the  mass  flow  is  reduced,  the  downstream 
expansions  disappear  first,  in  a  noticeable  step-wise  fashion  until  the  inlet 
expansion  ceases  to  be  present.  At  this  point,  it  is  found  that  violent  and 
unannounced  Impeller  surge  takes  place . 

From  the  theoretical  work  of  this  section  some  insight  can  be  gained  as 
to  the  nature  of  the  phenomena  taking  place  in  a  transonic  impeller  operated 
at  constant  speed  with  the  mass  flow  gradually  reduced.  Study  of  Figure  70 
suggests  that,  as  the  back  pressure  goes  up,  shock  MU  jumps  from  the 
oblique  to  the  normal  position  and  ail  points  R,  N,  T,  M.  etc.  become 
compressed  a  little  towards  the  upstream  direction.  With  further  pressure 
increase  shock  NT,  in  turn,  becomes  normal.  This  accounts  for  the  observed 
step-wise  change  in  pressure,  which  progresses  from  downstream  to  upstream. 
Finally,  the  Mach  number  at  R  is  no  longer  high  enough  to  permit  a  reflected 
oblique  shock  at  R,  and  shock  CR  itself  becomes  normal.  It  may  well  be 
that  the  presence  of  a  strong  normal  shock  at  the  inlet  causes  widespread 
boundary  layer  separation  there .  The  resulting  sudden  drop  in  impeller 
efficiency  would  bring  about  impeller  surge. 

it  is  clear  that  the  above  described  process  may  be  quite  efficient  up  to 
the  instant  when  boundary  layer  separation  occurs.  In  order  to  retard  this 
occurrence,  the  blade  design  muet  be  such  that  either  the  oblique  shock  pattern 
can  be  maintained  as  long  as  possible  or  efficient  operation  can  be  had  with 
a  normal  shock  at  the  inlet.  In  view  of  the  past  lack:  of  success  of  the  so-called 
3hock-in-rotor  type  of  compressor,  the  former  alternative  seems  to  be 
preferable . 
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SUMMARY 


The  basic  objective  of  the  present  work  is  to  formulate  basic  criteria 
for  the  design  of  centrifugal,  and  mixed-flow  impellers.  These  criteria 
must  be  general  ir  nature  and  must  be  so  selected  as  to  supply  design 
Information  in  the  specific  areas  where,  presently,  the  design  is  left  to 
judgement  and  to  the  experience  of  the  designer.  A  second  objective  is 
to  improve  the  understanding  of  the  internal  aerodynamics  so  as  to  advance 
the  state  of  the  art  of  high  performance  compressor  design. 

A  study  to  allow  selection  of  the  quantities  to  be  investigated,  indicated 
that  the  over -all  blade  loading,  tne  distribution  of  blade  loading,  and  the  exit 
cone  angle,  are  most  significant  quantities  in  that  their  effect  on  compressor 
performance  is  large  and  is  not  well  understood,  in  order  to  investigate 
these  effects  experimentally,  ten  impellers  were  designed.  Two  series  oi 
three  and  four  impellers,  respectively,  were  designed  to  investigate  the 
effect  of  impeller  length  on  performance  and  simultaneously  that  of  over-all 
blade  loading.  A  third  series  of  three  impellers  was  made  to  study  the 
effect  of  exit  cone  angle.  Within  each  series,  the  maximum  blade  loading 
of  each  impeller  Ls  varied  from  inducer  to  exit.  Each  impeller  was  so 
designed  as  to  isolate,  as  completely  as  possiole,  the  effect  of  the  primary 
variable  under  study.  The  characteristics  of  each  impeller  are  fully  described 
in  this  report  as  well  as  the  design  procedure  used  to  compute  these  character¬ 
istics  . 

Future  work  will  involve  fabrication  and  testing  of  the  three  most  promising 
impellers  as  an  initial  step  towards  attaining  the  desired  objectives. 
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Figure  1 


EFFECT  OF  IMPELLER  DECELERATION  RATIO  UPON 
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MEASURED  FLOW  COEFFICIENT  IN  TYPIC  \L  RADIAL  INLET  FLOW  PASSAGE  OF  COMPRESSOR 
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RELATIVE  MACH  NUMBER  ALONG  IMPELLER  LEADING  EDGE 
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ABSOLUTE  MACH  NUMBER 


tOBHPEHTIAl 


NOTE:  (INCLUDES  EFFECTS  OF  VORTICITY, 
VISCOUS  AND  BLADE  CLOGGING) 
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SHAPE  OF  LEADING  EDGES  -  THEORETICALLY  CALCULATED 
SHAPE  AND  FIRST  DESIGN  SHAPE 
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FRONT  VIEW  OF  IMFELLER 
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FRONT  VIEW  OF  IMPELLER 
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FRONT  MEW  OF  IMPELLER 
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TOTAL  LENGTH  FROM  IMPELLER  INLET  TO  DISCHARGE 


PRESSURE  DISTRIBUTION  ALONG  SHROUD  AND  HUB  FOR  IMPELLER  IIB 
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Figure  65 
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Inducer  Impe'lers  Figure  20  Figure  46 

Figure  47 

Lightly  loaded  Figure  13  Figure  49  Figure  52 

Inducer  Figure  21  Figure  50 

Figure  51 
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TABLE  H 

IMPELLER  DESIGN  CONDITIONS 


Inlet  Conditions 

oOF  and  29.92  In.  Hg 

Mass  Flow  W/e/5 

401  lb  per  minute 

Impeller  efficiency 

0.90 

impeller  Work  Input 

g 

73, 500  ft-ib  per  lb 

Slip  Factor 

0.81 

Flow  Factor 

See  Figure  4 

Relative  Deceleration  Ratio, 

wi 

0.375 

Tip  Speed 

1740  ft  per  sec 

RPM//0 

45,100 

Nominal  work  ratio  betweeiwihrcud 
and  hub,  (U  cul)s/(l4  cua  )h 

1.06 

BASIC  GEOMETRY  COMMON  TO  ALL  IMPELLERS 
Impeller  Mean  Exit  Diameter  8.83  inches 

Inducer  Tip  Diameter  to  Mean  Impeller  Exit  Diameter  Ratio  0.713 


COMMON  IMPELLER  OPERATING  CONDITIONS 


Total  Temperature  Rise  3 OOF 

Static  Temperature  Rise  206F 

Nominal  Maximum  Inlet  Relative  Mach  No.  1 .32 
Mean  Discharge  Absolute  Mach  No.  1.14 

Impeller  Total  Pressure  Ratio  6.11 

Impeller  Static  Pressure  Ratio  2.82 

Degree  ol  Reaction  0  535 
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Axial 

Coord. 

Z(in.) 

-5.800 

-5.400 

-5.000 

-4.600 

-4.200 

-4.000 

-3.800 

3.GGC 

-3.400 

-3.300 

-3.200 

-3.100 

-3.000 

-2.800 

-2.800 

-2.700 

-2.600 

-2.500 

-2.400 

-2.300 

-2.200 

-2.100 

-2.000 

-1.800 

-1.800 

-1.700 

-1.600 

-1.500 

-1.400 

-1.300 

-1.200 

-1.100 

-1.000 

-0.900 

-0.800 

WADC  TR56 


TABLE  m 

INLET  COORDINATES 


R;idiu3 

Radios 

to  Hub 

to  Shroud 

Rh(in.) 

«.(*»-> 

4.062 

3.750 

3.418 

3.086 

2.754 

2.588 

2.422 

2.256 

2.090 

2.008 

4 

1.926 

1.848 

1.776 

1.708 

1 .645 

1.588 

1.534 

1.482 

1.435 

1.390 

1.348 

1.313 

1.280 

1.249 

1.221 

1.194 

1.170 

1.149 

1.132 

4.360 

1.116 

0.634  (radius) 

1.105 

3.451 

«  OA  rr 

1  VT»  t 

3.394 

1.090 

3.344 

1.084 

3.301 

1.060 

3.261 
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TABLE  m  (Coni’ d) 


Axial 

Radius 

Radius 

Coord. 

to  huu 

to  tihroud 

Z(in.) 

Rh(la.) 

-Vto> 

-0.700 

1.078 

3.229 

-0.600 

1.078 

2.200 

-0.500 

1.075 

3.183 

*0.400 

1.075 

3.168 

-0.300 

1.075 

3.157 

-0.20C 

1.075 

3.152 

-0.100 

1.075 

3.150 

-0.000 

1.075 

3.150 
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TABLE  IV 

IMPELLER  SERIES  I  -  HUB  AND  SHROUD  COORDINATES 


Axial 

Radius 

Radius 

Radius  to 

Coord . 

to  Hub 

to  Shroud 

Blade  Tip 

Z(in.) 

Rh(ln.) 

«s(in.) 

-0.000 

1.075 

3.150 

3.135 

0.100 

1.080 

3.150 

3.135 

0.ZUQ 

1.391 

3.150 

3.135 

0.300 

1 . 103 

3.150 

3.135 

0.400 

1.119 

3.150 

3.135 

0.500 

1.136 

3.151 

3.136 

0.600 

1.158 

3.152 

3.137 

0.700 

1.162 

3.156 

3.141 

0.800 

1.215 

3.165 

3.150 

0.900 

1.251 

3.178 

3.162 

1.000 

1.290 

3.192 

3.176 

1.100 

1.335 

3.213 

3 .197 

1.200 

1.384 

3.237 

3.220 

1.300 

1.438 

3.265 

3.248 

1.400 

1.478 

3.297 

3.280 

1.500 

1.561 

3.334 

3.317 

1.600 

1.634 

3.3.75 

3.358 

1.700 

1.714 

3.420 

3.403 

1 .300 

1.804 

3.472 

3.454 

1 .900 

1.906 

3.528 

3.510 

2.900 

2.020 

3.589 

3.571 

2.100 

2.148 

3.657 

?  638 

2.200 

2.2S4 

3.732 

3.712 

2.300 

2.446 

3.816 

3.796 

2.100 

2.611 

3.911 

3.889 

2.500 

2.782 

4.020 

3.998 

2.600 

2.954 

4.128 

4.104 

2.700 

3.125 

4.250 

4.226 
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TABLE  IV  (Cont'd) 


Axial 

Radius 

Radius 

Radius  to 

Coord. 

to  Hud 

to  Shroud 

Blade  Tip 

Z'in.) 

R.  (in./ 

« 

Rg(i*.> 

RjUn.) 

2.800 

3.209 

4.378 

4.352 

2.000 

3.468 

4.517 

4.489 

2.930 

4.560 

2.943 

4.548 

3.000 

3.846 

3.  ICO 

3.819 

3.200 

4.000 

3.300 

4.166 

3.360 

4.270 
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TABLE  V 

IMPELLER  SERIES  II  -  HUB  AND  SHROUD  COORDINATES 


Axial 

Radius 

Radius 

Radius  to 

Coord. 

to  Hub 

to  Shroud 

Blade  Tip 

Z(in.) 

Rh(in.) 

Rs(in.) 

-1.2C0 

1.105 

3.451 

-1.100 

1.097 

3.394 

-1.000 

1.090 

3.344 

-0.900 

1.087 

3.298 

-0.800 

1.084 

3.257 

-0.700 

1.083 

3.223 

-0.600 

1.083 

3.195 

-0.500 

1.083 

3.169 

-0.400 

1.086 

3.148 

*0.300 

1.090 

3.131 

-0.200 

i.096 

3.117 

-0.100 

1.103 

3.106 

0.G00 

1.112 

3.096 

0.100 

1.120 

3.089 

0.200 

1.130 

3.083 

0.300 

1.139 

3.080 

0.400 

1.149 

3.07? 

0.500 

1.160 

3 .075 

0.600 

1.171 

3.075 

0.700 

1.184 

3.076 

3.061 

0.800 

1.199 

3.078 

3.063 

0.900 

1.215 

3.083 

3.067 

1.000 

1.234 

3.091 

3.075 

1.100 

1.256 

3.102 

3.086 

1.200 

1.282 

3.119 

3.103 

1.300 

1.315 

3.141 

3.124 

1.400 

1.352 

3.167 

3.150 

k  «  vwu 

1.395 

3.198 

3.181 

1.600 

1.440 

3.235 

3.218 

1.700 

1.491 

3.283 

3.264 

1.800 

1.542 

3.337 

3.318 

1.900 

1.597 

3.400 

3.381 
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TABLE  V  (Cont’d) 


Axial 

Radius 

Radius 

Radius  to 

Ccord. 

to  Hub 

to  Shroud 

Blade  Tip 

Z  (in.) 

Rh(in.) 

Rs(in.) 

Rx(ln.) 

2.000 

1.656 

3.473 

3.453 

2.100 

1.718 

3.553 

3.531 

2.200 

1.782 

3.646 

*  «24 

2.300 

1.857 

3.749 

d.iiiS 

2.400 

1.940 

3.863 

3.839 

2.500 

2.044 

3.984 

3.958 

2.600 

2.175 

4.115 

4.085 

2.700 

2.345 

4.245 

4.219 

2.800 

2.608 

4.378 

4.352 

2.900 

3.000 

4.516 

4.487 

2.930 

4.560 

2.943 

4.548 

3.000 

3.357 

3.100 

3.680 

3.200 

3.959 

3.300 

4.150 

3.360 

4.270 
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TABLE  Vi 

IMPELLER  I  —  BLADE  COORDINATES 


Axial 

Coozd. 


ZQ°.) 


Angular 


Coord, 

Xdeg) 


_ Elide  Tlilck ztesztB 

Shroud 

Nofrm.  Ta  rig.  Norm. 


3uF - 

Tang. 


-0.700 

-23.040 

0.0280 

-0.650 

-20.970 

0.0280 

-0.600 

-18.920 

0.0300 

-0.550 

-18,880 

G.GS60 

-0.5000 

-14.870 

0.0440 

-0.450 

-12.926 

0.0475 

-0.400 

-11.330 

0.0200 

0.0368 

0.0505 

-0.350 

-  e.671 

0.0230 

0.0421 

0,0540 

-0.300 

-  8.160 

0.0255 

0.0464 

0.0570 

-0.250 

-  6.707 

0.0270 

0.0486 

0.0600 

-0.200 

-  5.285 

0.0295 

0.C525 

0.0625 

-0.150 

-  3.908 

0.0320 

0 . 055£ 

0.0650 

-0.10C 

-  2.570 

0.0345 

G.C6GC 

0.0675 

-0.050 

-  1.273 

0.0365 

0.0627 

0.0700 

0.090 

0.000 

0.0385 

0.0653 

0.0720 

0.050 

1.244 

0.0405 

0.0668 

0.0745 

0.100 

2.460 

0.0425 

0.0701 

0.0780 

0.150 

3.659 

0.0445 

0.0722 

0.0785 

0.200 

4.830 

0.0465 

0 . 0738 

0.0805 

0.250 

5.982 

0.0485 

0.0758 

0.0825 

0.300 

7.105 

0.0505 

0.0788 

0.0845 

0.350 

3.190 

0.0525 

0.079° 

0.0860 

0.400 

9.250 

0.0545 

0.0823 

0.0880 

0.450 

10.286 

0.0665 

0.0840 

0.0900 

0.500 

11.300 

0.0580 

0.0855 

0.0915 

0.550 

12.294 

0.0590 

0. OP  fid 

0.0930 

0.600 

13.260 

0.0595 

0.0845 

0.G945 

0.650 

14.192 

0.0600 

0.0835 

0.0965 

0.700 

15.100 

0.0600 

0.0825 

0.0960 

0.750 

15.971 

0.0810 

0.0990 

0.800 

16.820 

0.0800 

0,1005 

0.850 

17.660 

0.0785 

0.1013 

0.900 

18.430 

0.0770 

0.1030 

0.950 

19.175 

0.0760 

0.1040 

1  .GOO 

19.880 

0.0600 

0.0750 

0.1050 
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0.0373 
0.0373 
0.0400 
0.0480 
0.0530 
0.0570 
0.0300 
0.0630 
0.0055 
0.0S85 
0.0710 
0.0735 
0.0750 
0.0770 
0.0795 
0.0815 
0.0335 
0.0855 
0.0880 
0.0900 
0.0915 
0.0935 
0.0950 
0.0970 
0 .0930 
0.0995 
0.1005 
0.1020 
0.1030 
0.1040 
0.1050 
0.1060 
0.1070 
0.1080 
0.1090 
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TABLE  VI  (Cont’d) 


Axial 

Angular 

Blade  T'lickneft&es 

Coord. 

Coord. 

Siu  cut* 

Hub 

Z  (In.) 

tldej£)_ 

Norm. 

Tang. 

Norm. 

Tang. 

1.100 

21  180 

0.0600 

0.0725 

0.1070 

0.1100 

1.200 

22.350 

| 

1 

0.0700 

0.1080 

0.1105 

1.300 

23.405 

0.0680 

0.1085 

0.1105 

1.400 

24.360 

0.0665 

0.1080 

0.1095 

1.500 

25.220 

0.0650 

0.1060 

0.1075 

1 .600 

25.980 

0.0640 

0.1025 

0.1045 

1.700 

26.640 

0.0635 

0.0990 

0.1005 

1.800 

27.210 

0.0630 

0.0960 

0.0970 

1.800 

27.700 

0.0620 

0.0935 

0.0940 

2.000 

28.145 

0.0620 

0.0920 

0.0920 

2.100 

28.585 

0.0615 

0.0910 

0.0910 

2.200 

29.025 

0.0610 

0.0900 

0.0900 

2.300 

29.400 

0.0810 

0.0900 

0.09C0 

2.400 

29.760 

0.0610 

0.0900 

0.0900 

2.500 

30.090 

0.0605 

0.0900 

0.0900 

2.600 

30.390 

0.0600 

0.0900 

0.0900 

2.700 

30.670 

0.0585 

0.0590 

0.0900 

0.0900 

2.800 

30.940 

0.0560 

0.0565 

0.0900 

0.0900 

2.900 

31.195 

0.0515 

0.0515 

0.0900 

0.0900 

2.943 

0.0480 

0.0480 

3.000 

31.470 

0.0499 

0.0499 

0.0880 

0.0880 

3.100 

31.730 

0.0534 

0.0534 

0.0850 

0.0850 

3.200 

31.365 

0.0569 

0.0569 

0.0800 

0.0850 

3.300 

32.175 

0.0604 

0.0604 

0.0710 

0.0710 

3.360 

32.300 

Q.0625 

0.Q625 

0.0630 

0.0630 
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TABLE  VJI 

IMPELLER  LB  ---  BLADE  COORDINATES 


Axial 

Angular 

Blade  Yhic 

knessej 

Coord. 

Coord. 

Shroud 

Hub 

Z(in.) 

*(deg) 

Norm. 

Tang. 

Norm . 

Tang. 

-0.700 

-27.9138 

0.0280 

0.03G1 

-0.650 

-25.663 

0.0280 

0.0361 

-0.600 

-23.4376 

0.0300 

0.0386 

-0.550 

-21.2376 

0.0380 

0.0487 

-0.500 

-19.0630 

0.0440 

0.0563 

-0.450 

-16.9029 

O.o475 

0.0605 

-0.400 

-14.7946 

0.0505 

0.0641 

-0.550 

-12.7351 

0.0540 

0.0683 

-0.300 

-10.7313 

0  0570 

0.0713 

-0.250 

-  8.7856 

0.0600 

0.0744 

-0.2QQ 

-  6.9004 

0.0625 

0.0763 

-0.150 

-  5.0828 

0.0650 

0.0783 

-0.100 

-  3.330* 

0.0675 

0.0804 

-0.050 

-  1  6349 

0.0700 

0.0825 

0.000 

0.0000 

0.0200 

0.0410 

0.0720 

0.0838 

0.050 

1.5673 

0.0245 

0.0490 

0.0745 

0.0860 

0.100 

3.0847 

0.0285 

0.0569 

0 . 0760 

0.0870 

0.150 

4.5588 

0.0330 

0.0648 

0.0785 

0.0892 

0.200 

5.9922 

0.0370 

0.0715 

0.0805 

0.0906 

0.250 

7.3868 

0.0415 

0.OY80 

0.0825 

0.0926 

0.300 

8.7451 

0  0455 

0.0831 

0.0845 

0.0946 

0.350 

10.0691 

0.0500 

0.0890 

0.0860 

0.0959 

0.400 

11.3502 

0.0536 

0.0931 

0.0880 

0.0979 

0.450 

12.6009 

0.0565 

0.0359 

0.0900 

rs  poo  n 

0.500 

13.8126 

0.0580 

0.0964 

0.0915 

0.1012 

0.550 

14.9926 

0.0590 

0.0962 

0.0930 

0.1025 

0.600 

16.1420 

0.0595 

0.0957 

0.0945 

0.1040 

0.650 

17.2627 

0.0600 

0.0950 

0.0965 

0.1059 

0.700 

18.3565 

i 

0.0934 

0.0980 

0.1071 

0.750 

l''. 4247 

0,0924 

0,0990 

0.1080 

0.800 

20.4634 

0.0912 

0.1005 

0.1091 

0.850 

21.4845 

0.0897 

0.1015 

0.1099 

0.900 

22.4943 

0.0884 

0.1030 

0.1112 

0.950 

23.4889 

0.0875 

0.1040 

0.1120 

1 .000 

24.4593 

0.0600 

0.0853 

0.1050 

0.1129 
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TABLE  VH  (Cont’d) 


Axial 

Angular 

Blade  Ttiicknesses 

Coord. 

Coord. 

Shroud 

~Hub 

Z(ln.) 

Norm. 

Tang. 

Norm. 

Tang. 

1.100 

26.3367 

0.0600 

0.0839 

0.1070 

0.1144 

1.200 

28.1197 

0.0818 

0.1080 

0.1150 

1.300 

29.8185 

0.0797 

0.1085 

0.1152 

1.400 

31 .4391 

0.0776 

0.1080 

0.1144 

1.500 

32.9869 

0.0760 

0.1060 

0.1121 

1.600 

34.4760 

0.0741 

0.1025 

0.1081 

1.700 

3L-  3002 

0.0724 

0.0990 

A  .  A  .  A 

U.iUlV 

1.800 

37.2585 

i 

i 

0.0709 

0.0© 60 

0.1002 

1.000 

38.5461 

1 

0 . 0694 

0.0S35 

0.0973 

2.000 

39.7540 

0.0600 

0.0680 

0.0920 

0.0952 
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TABLE  Vin 


IMPELLER  n  —  BLADE  COORDINATES 


Axial 

Angular 

Blade  Thicknesses 

Coord. 

Coord. 

Shroud 

Huu 

Z(ln.) 

0  (d€<n*) 

Norm. 

Tang. 

Norm. 

Tang. 

0.1G0 

0.000 

0.0280 

0.0371 

0  200 

4.134 

0.0305 

0.0383 

0.250 

6.031 

0.0343 

0.0428 

0.300 

7.903 

0.036? 

0.0455 

0.350 

9.739 

0.0362 

0.0470 

0.400 

11.532 

0.0304 

Q .0483 

0.450 

13.278 

0.0407 

C .0496 

0.500 

14.975 

0.0420 

0 .0507 

0.550 

16.622 

0.0206 

0.0412 

0.0431 

0.0517 

0.600 

18.209 

0.0202 

0.0395 

0.0442 

0.0525 

0.650 

19.721 

C.0201 

G.G382 

0.04:54 

0.0532 

0.700 

21.155 

0.0203 

0.0376 

0.04167 

0.0539 

0.750 

22.499 

0.0214 

C.CG3I 

0.0479 

0.0545 

0.800 

23.774 

0.0230 

0.0392 

0.0(90 

0.0552 

0.850 

24.984 

0.0248 

0.0410 

Q.0502 

0.0559 

3  9CC 

26.124 

0.0270 

0.0425 

0.0513 

0.0565 

0.95C 

27.204 

0.0291 

0.0442 

0.0525 

0.0574 

1.000 

28.224 

0  0311 

0.0457 

0.0538 

0.0584 

1.100 

30.134 

0.0346 

0.0480 

0.0560 

0.0601 

1.200 

31.868 

0.0374 

0.0495 

0.0584 

0.0620 

1.300 

33.457 

0.0393 

0.0500 

0.0607 

0.0640 

1.400 

34.922 

0.0407 

0.0502 

0.0628 

0.G658 

1.500 

36.264 

0.0420 

0.0505 

0.0641 

0.0675 

1.600 

37.494 

0.0435 

0.0510 

0.0668 

0.0693 

1.700 

38.622 

0.0450 

00517 

0.0687 

0.0710 

1.800 

39.649 

0.0466 

0.0526 

0.0704 

0.0726 

1.900 

40.575 

0.0483 

0.0537 

0.0722 

0.0743 

2.000 

41.414 

0.0500 

0.0550 

0.0744 

0.0764 

2.100 

42.164 

0.0517 

0.0561 

0.0780 

0.0798 

2.200 

42.837 

0.0534 

0.0572 

0.C852 

0.0848 

2.300 

43.449 

0.0552 

0.0586 

0.0907 

0.0921 

2.400 

44.004 

0.0589 

0.0598 

0.1000 

0.1012 
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TABLE  VIE  (Cont'd) 
IMPELLER  n  —  BLADE  COORDINATES 


Axial 

Angular 

Blade  Thicknesses 

Coord. 

Coord. 

Shroud 

Hub 

am.) 

ti*Zg£l 

Norm. 

Tang. 

Norm. 

2.500 

44.489 

0.0584 

0.0607 

0.1121 

0.1131 

2.600 

44.944 

0  0597 

0.0615 

0.1251 

0.1260 

2.700 

45.339 

0.0597 

0.0610 

0.1328 

0.1332 

2.800 

45.674 

0.0577 

0.0585 

0.1356 

0.1330 

2.900 

45.924 

0.0485 

0.0488 

0.1335 

C\  me 

3.000 

46.048 

0.0463 

0.0455 

0.1286 

0.1286 

3.100 

46.164 

0.0500 

0.0502 

0.1203 

0.1203 

3.200 

46.312 

0.05r  j 

0.0538 

0.10 to 

0.1061 

3.300 

46.504 

0.05  2 

0.0575 

0.0838 

0.0840 

3.360 

46.661 

0.05.5 

0.0598 

0.0595 

0.059S 
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TABLE  IX 

LEADING  EDGES  FOR  IMPELLERS  I,  IB,  D 


Shape  cf  initial  leading  edge  profile 


Axta'j 

Radius  to  point  on  leading  edge 

Cccrdins  tee 

Z  in. 

Impeller  I 

Impeller  IB 

Impeller  n 

-0.70 

1.075 

1 .075 

-0.65 

1.740 

1.235 

-0.60 

2.400 

1.400 

-0.55 

«.  aid 

i  .570 

-C.5Q 

2.785 

1.750 

-0.45 

2.95 

1 .895 

-0.40 

3.135 

2.060 

-0.35 

2.230 

-0.30 

2.400 

-0.25 

2.520 

-0.20 

2.540 

-0.15 

2.765 

. -0.10 

2.885 

-0.05 

3.01 

-0.0 

3.135 

1.120 

+0.05 

1.600 

+0.10 

2.200 

+0.15 

2.205 

+0.20 

.2.390 

+0.25 

2.490 

+0.30 

2.580 

+0.35 

2.675 

+0.40 

2.770 

+0.45 

2.870 

*•0.50 

2.665 

+0.55 

3.060 
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NOMENCLATURE 


Sonic  verity  at  temperature  -  T 
Area 

Radius  of  sphere  and  cylinder 

Absolute  velocity  of  fluid 

Specific  heat  at  constant  pressure 

Inducer  diameter 

Impeller  moan  exit  diaxuater 

Blade  force  upon  fluid,  in  radial  direction 

1  “  ( -71k—)  * 

°t 

Gravitational  constant 
Change  in  total  head 
Enthalpy 

Mechanical  equivalent  of  heat 
Ratio  of  specific  heats 


m,  a 


cmjj  aver./^mi,  s 

Rotative  speed,  revolutions  per  sec 

System  of  coordinates  relative  to  the  streamlines  in  the 
meridional  plane;  mis  along  the  streamline,  increasing  from 
Impeller  inlet  to  *xifc  n  is  normal  to  the  streamline,  increasing 
from  shroud  to  hub 
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ft 
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P 
Q 
r 
R 
R 
8 
T 
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u 
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U 
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W 
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z 

Z 
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_ L  _ 4 

_ _  _  . -  —  


Polytropic  exponent 


i 

Impeller  volume  flow  under  inlet  conditions 
Radial  coordinate  measured  from  axis  of  rotation 
Perfect  gas  constant  for  air 

» 

Streamline  radius  of  curvature  in  the  meridional  plane 
Entropy 

i 

Temperature  R 

* 

Tangential  blade  thickness 
iuipeller  tangential  velocity  -  ro> 

Impeller  mean  exit  tip  speed 

Uniform  flow  at  -  •»  for  a  cylinder  or  a  sphere 

Impeller  relative  velocity 

Compressor  mass  rate  of  flow 

Axial  coordinate,  increasing  from  the  impeller  inlet  to  the 
impeLer  exit 

Number  of  blades 
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(  )  aver 

(  )  B 

(  )C 
(  )  H 

(  )I 
(  )  i® 

(  ) 

(  )  » 

(  )P 
(  )  Pr 

(  )  r 
(  >♦ 

(  )  8 
i  )  suet 
(  )u 
(  )  x 
(  ) 2 
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SUBSCRIPTS 


Average 

Blade 

Compressor 

Hub 

Impeller 

Isentropic 

Refers  to  the  meridional  component 

Refers  to  a  component  normal  to  the  streamline  in  the 
meridional  plane 

Potential  flow 

Applies  to  the  blade  pressure  surface 
Radial  component 
Along  a  streamline 
Shroud 

Supplies  to  the  blade  suction  surface 
Tangential  component 
Impeller  blade  tip 
Axial  component 
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3 

5 

n 

€ 

A 

V 

I 

4 

t 

p 

a 

* 

!t> 


GREEK  SYMBOLS 

Angie  between  the  axis  of  rotation  and  .in  tangent  to  a  stream 
line  in  th€  meridional  plane 

Angle  between  the  relative  flow  direction  aod  the  meridional 
plane 

Pressure/ 14. 7 
Efficiency 

Viscous  clogging  coefficient 
Temperature  coefficient 
Specific  Speed 
Flow  parameter 

Angular  coordinate  in  a  plane  normal  to  the  impeller  axis 
of  rotation 

Stream  function,  such  that 

a*  p 

— *-  =  r  —  w. 

dr  P„  * 

Mass  density 

Geometrical  parameter  -fy~  a 

Temperature/ 520,  Polar  coordinate  for  cylinder  and  sphere 

Rotational  speed,  rad  per  sec 

Vorticlty 
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APPENDIX  A 

DERIVATION  OF  RELATIONS  CONTROLLING 
TEE  IMPELLER  EXTERNAL  OVER- ALL  GEOMETRY 


(1)  IMPELLER  INIJ5T  TO  EXIT  DIAMETER  RATIO 

(«)  flalfttlca  between  specific  speed  a ad  iialet  to  wdt  diameter 
ratio 

t» •  N^/Q 

(gAH)% 

Allowing  N  *  jk 

gAH*  i^AU* 

(cm,)ovar*  K  c^„ 

^^iirr*40”^ 

There  follows 

u  =  0.282yRT  I 

In  this  equation: 

(1)  tan  can  be  evaluated  for  the  condition  of  mini¬ 
mum  inducer  tip  relative  Mach  Number 
(See  part  B  below. ) 


(2)  For  conventional  impellers  the  ratio  /j{f 
ia  never  very  far  from  1 . 0.  <V>*A 


(3)  Impeller  specific  speed  is  then  directly  related  to 
the  ratio  of  impeller  inlet  to  exit  diameters. 
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(b)  Value  of  &$•  for  minimum  Laid  relaLve  velocity 

This  condition  if?  ve»-r  nearly  that  of  minimum  mdncer 
?:p  relative  Mach  Number  It  will  be  worsted  out  for  coostant 
rpm  and  constant  mass  rate  of  fkrj 


(Mv1^  ds*  ♦ 


{ Wc/g  j*  )* 
(K  f  x  }* 


Introducing  A*  =  tt 


Rq gating  the  right  fcand  side  to  aero  and  solving  for  d5 


ir( a* — )*» . a  cap*' ti *  1  % 

°s  SjfKf5/k  l  /*  A*  d  $  A/A*)  -5 


This  relation  may  now  be  introduced  into  an  expression  for  *an  n_ 


2{ 


.  A 

1  A*  di  A/A*rJ 


(A-2) 


The  ritijt-hand  side  is  a  function  of  the  average  inlet  Mach 
Number  only.  It  is  so  plotted  on  figure  2- 
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Expression  for  specific  speed  in  terras  of  rareentlonal 
compressor  coordinates 


It  Is  convenient  to  express  dimensionless  specific  speed 
in  terms  of  more  conventional  compressor  coordinates 


N^/5 

y  S 

(gAH)% 

N  = 

RPM/60 

Wc. 

Q  * 

6°5A>0( 

where 

wc 

where 


0.2185 


v*feas, 

801 


la  in  Ibs/min. 


gAH  s  JgCp  (AT0)j# 
*  6060 


3  150  000  ®oi  [(£tt)  is* 

pw 


k-l 


Combining  all  this 


J _  RPM  (Wcv^ 

*  =  0.9I3XI07  Je^  '  So. 


2095 

10* 


K 


0*8m 
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(d)  Expression  for  the  ratio  of  inducer  tip  diameter 
to  unpeuer  tip  diameter. 


k-l 


In  equation  A- 5,  we  may  replace  f(^”“)k  - 1 1 


199  5.15  x  10* 


\rj 


Li* 

I  6o. 


Equation  A- 1  and  A- 5  in  this  form. 


0.282 


/k7v  \ 

(ijji)  M  v 'tan  3s 


2.095  f  J _ l_ 

*0*  Vfc  0 


( 


*/T _ 

So,  .  i 

3.15  x  10* 


! 


I  dot 


and  after  simplification. 


r_J_  Kf  ~  I 
'v'So.  ./ton  3* 


So* 


Jinches  = 


( A-6) 


This  equation  conveys  essentially  the  same  infor 
roation  as  equation  A- 1 .  Its  form  is,  however,  more 
convenient  and  it  is  plotted  on  Figure  3 . 
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(?)  IMPELLER  EXIT  WIDTH 

(a)  State  of  fluid  at  impeller  exit.  The  state  of  the  Held 
anywhere  in  the  impeLer  crx  be  expressed  in  term*  t><  the 
impeller  tangential  velocity  and  tne  fluid  relative  velocity 
from  the  first  law. 


Trors  the  geometry 


wa  K  c*  +  u*  —  2ucu 


Combining 


*  w«  -  u6 


x 

Tw 


H  /  u  * 

T{cl 


(A-7) 


which  is  useful  in  calculating  impeller  temperature*.  Im¬ 
peller  pressures  end  densities  ars  computed  by  assuming 
isentrcpic  relations  to  bold  in  the  inlet,  upstream  of  the  lead' 
lag  edge  and  adiabatic  flow  with  losses  to  exist  within  the 
impeller.  On  this  the  density  ratio  follows  as 


*<*«%.)<£ 

■  (Aov/^) [•  • ♦  x  < Sg  -  5r ' ' 'Vl  <A' 85 
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Similarly  for  pressure 


Ri 


PM 


Pot 


T.ovl  2  ‘o' 


Of 


-Sr>]} 


%H 


(A-9) 


la  this  equation 


To. 
Tr  ov 


Pw_ 

P.OV 


Pot 

P*cn 


are  functions  o i  the  average  inlet  Mach  Number.  The  quantity 
o  is  a  polytropic  exponent,  which  differs  from  k  in  that 
it  mahes  allowance  for  estimated  impeller  losses.  It  is  related 
to  the  impeller  efficiency  in  J*e  following  maimer. 


314 

P* 

II 

-i) 

but 

(TotVs  (Ta)is  T.  , 

To.  ~  T.  To.  1 

Tc  *  . 

tT1* 

]o»_  .  _Tjl  T« 

To* 

To.  T  i  To. 

To. 

Here  (T“)js  -  » 

!  t 

and 


Pt 

P, 


(PtfvS 

P, 


combining 


T. 

To. 


^H(pT)K  -<p7*  "  1  1  ‘  ’i  (a-io) 


K-i 


n-1 


For  design  purpose,  the  value  of  n  may  be  computed  from 
this  equation  for  desired  values  of  wheel  reaction  end  wheel 
inlet  velocities,  and  for  estimated  impeller  efficiencies. 
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(b)  Impeller  Viscous  flow  coefficient 

(1)  Original  definition 

In  the  design  procedure  oi  this  report,  the  viocoue 
flow  coefficient  may  be  defined  as  the  ratio  of  the  impeller 
exit  width  to  the  theoretical  width  that  would  satisfy  mass 
flow  impeller  efficient  and  the  exit  velocity  triangle.  The 
velocity  triangle  takes  into  consideration  all  the  im¬ 
peller  losses  without  separating  impeller  win&ge  and  internal 
losses  which  thus  affect  both  the  leaving  total  pressure  and 
the  leaving  tangential  velocity  component.  Thus  the  viscous 
flow  coefficient  is  more  or  less  a  catch- ail  term  seeking  to 
accost  all  at  cn.ce  for  the  impeller  external  losses  and 
boundary  layer  development  as  well. 

Recent  experimental  work  at  AiResearch  shows  that 
this  coefficient,  as  defined  above,  can  be  either  smaller 
or  greater  than  1.0.  That  this  is  possible  ip  made  clear 
if  one  realize;}  that  consistency  must  ho  had  between  the 
impeller  efficiency  arid  the  flow  factor.  A  flow  factor  is 
necessary  only  insofar  as  it  relates  the  energy  profile  and 
the  mass  flaw  profile  at  the  exit.  This  factor  is  necessarily 
close  to  unity,  except  for  grossly  distorted  velocity  distri¬ 
butions.  The  introduction  of  the  windage  loss  into  this 
definition  can  only  be  justified  from  the  standpoint  of 
simplification  of  the  design  procedure. 

(2)  Exact  definition 

Two  rigorous  definitions  of  an  exit  flow  coefficient 
can  be  formulated.  Only  the  second  one  however  is  really 
useful  for  design  purposes. 
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The  exit  flow  coefficient  may  be  defined  as  the 
ratio  of  the  actual  to  the  theoretical  flow.  Theoretical 
flow  nay  be  visualized  as  that  calculated  from  a 
theoretical  lsentroptc  impeller  of  the  same  aerodynanic 
work  1  input  and  the  same  degree  of  reaction.  Such  a 
definition  is  however  quite  trivial.  Lillie  usefulness 
can  be  made  of  tne  restating  fkm  coefficient  which  farther 
has  the  disadvantage  of  having  the  rather  lot/  numerical 
value  of  0.  6-  0.  7. 

A  useful  definition  at  flow  coefficient  is  that  which 
is  prescribed  by  the  necessary  consistency  which  must 
exist  between  flow  and  energy  if  an  averaged  two-dimension¬ 
al  profile,  as  is  present  at  the  impeller  discharge,  is  to 
be  represented  by  a  single  velocity  triangle.  In  simple 
language,  it  is  a  factor  required  by  the  'act  that  the 
meridional  component  of  the  integrated  exit  absolute 
velocity  is  different  from  the  integrated  value  ol  the 
meridional  component  of  this  velocity.  TLe  resulting  flow 
coefficient  is  obviously  a  function  cf  the  shape  of  the 
velocity  profile.  For  a  good  exit  profile  at  the  discharge 
of  an  impeller  of  this  type,  it  is  about  0. 96.  This  co¬ 
efficient  is  the  only  one  ihat  needs  to  be  taken  into  account 
in  fixing  the  impeller  discharge  width. 
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APPENDIX  B 

EQUATIONS  FOR  THE  AXI-SYMHJiTRIC 
FLOW  IF  IMPELLER  PA3SAGES  - 
EQUATIONS  FOR  DESIGN  PROCEDURE 


Below  are  given  the  equations  ccntrolling  the  motion  o i  she  fluid 
Inside  the  impeller.  They  are  deriveiJ!  for  the  axi- symmetric  ease 
corresponding  to  »n  infinite  number  or  blades,  and  for  isrritcid,  in¬ 
compressible  fluid. 

I  GENERAL  CASE 

Derivation  of  the  particular  form  at  the  general  equation  used  in 
the  d©3ign  procedure  will  be  made  using  as  a  starting  point  the  vorticity 
component  in  the  meridional  plane. 

dwm  +  _  _  tan ft  rd(f Cu)  1  iB-t) 

dn  R  r  cos  a  *  dr  ■*  Z  ~  constant 


"be  vector  representing  the  total  vorticity  is  a  vector  directed 
along  the  blade  and  along  a  line  erf  constant  angular  momentum  in  the 
meridional  plane.  The  axial  com  xxient  of  this  vorticity  is 

J_  Mrcu) 

V  r  r  dr 

which  is  nothing  more  than  the  definition  of  vorticity  in  radial  coordinates 
(ref.  12).  The  tangential  component  is 


_  tanft  r  d(rcu)  *| 
r  cos o  l-  dr. 

Equation  (B-l)  may  also  be  derived  directly  from  the  Lorenz  equation*. 


Now  returning  to  equation  (El- 1) 

rcu  =  r*«  -  rw~4  ion  Q 

„  I  dp 

=  ~  r*m  77 

Z  -  cdiiiiQnt 


r  dtfcu)i 

l7r~J 


-  r  tan  /3 


~d  r 
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The  condition  of  radial  blade  element*  Is  Introduced  In  tbe  form 


rton£] 

LrcosoJ 


constant 


Z  *  constant 


7^-  =  -  tapi  ft  cos*  9  tan  a  ^7- 

0i  a» 

La  addition 


COS*/3  .  0 

+  ~  - —  ten  ft 


a"m 

sin  a  - 

dr 

dm 

dn 

d  a 

do 

sin  a  - 

do 

dr 

dm 

dn 

COSO 


CO&  o 


(B-3) 


(B-4) 


(B-5) 


Combining 


2  i  *  f] 


CB-'l,  '2,  -3,  -4,  -5) 


d  r  R 


ton  ft  r 
rcoso  l 


,  de  da  .  cos /3  ~ I 

( - —  sin  a  -  - —  cssa)  +  — —  sin  ft  I  - 

dm  dn  r  J 

«]-•»  tan  fl} 


dw^ni 


|  2  r  1#  ~ s,n  ^  005  ^  ton  a 

nft  .  *.1  .  a 

-j—  sin  pj-  r  ton  p  t“T7“  sin  c 


sin  e 


which  simplifies  to 
a*m 


—  (I  ♦  ton*  ft)  +  wm  [ tan*  3  tan*  a)-  tan*  £  ton  a 
n  1  k 


ztortei  «  0 

—  -  =  ten  ft  tan  a 

rcosa  J 


<*>  _  ton/S 

d  m  cos  a 
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H  SPECIFIC  LOCAL  REGIONS 

In  specific  regions  of  the  Impeller,  the  general  equation  can  be 
expressed  in  simplified  forms  suitable  for  tend  computation. 

(a)  Inducer  portion  of  impeller. 

Starting  with  equation  B-l  and  introducing  the  condition  of  radial 
blade  elements  in  terms  of  the  shrewd  conditions 


r±'£a>l 

L  dr  J 


dr  -*  Z  =  constant 


(B-lOO) 


in  the  Inducer 


COSO 

COSO, 


C2«  » 

,an  a=  -sn; 

3  1.0 


4-  *  2r«-  -f-“  [r1  ~r~  *  2rwJ 

or  rt  L  dr 


Eq  (fa- 100)  becomes 


ton 
cos  os 


{2-n 


-  tan  +  2wm  -i  ]} 

As  a  further  simplifying  assumption, 

£  wrn  a  jym 

dr  dn 

in  the  inducer  section  and  the  final  equation  becomes 

_  2  tQn*&  Xl  s  -  pm  ton  ft  JL 

coso,  r,*J  cos  a,  r« 


ton  §% 

cosaK  r* 


(B- 102) 
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(b)  Downstream  porting  of  impeller 

» 

In  tbs  downstream  portions  of  the  impeller,  simplifying 
assumptions  other  than  those  in  the  inducer  portion  may  be  made 
for  easier  handling  of  the  equation.  In  view  cf  the  fact  that  little 
distortion  of  the  meridional  streamlines  occurs  in  the  ao*a~streara 
portion,  lit  is  pertinent  to  assume  that  their  radii  of  curvature  will 
remain  nearly  equal  to  the  radii  of  curvature  of  the  potential  flow 
through  the  corresponding  bladeless  passage.  Starting  from 
equation  (B-i01) 


i!to  +  J9n 

in  R 


c*»  ®s 


[*"  T. 


a  i  r *  ^wr 

-  fan  0,  < ,-T  77- 


?7  *] 


(8-IOI) 


Here,  ws  express  the  meridional  velocity  in  terms  of  meridional 
velocity  for  sore  vorticiiy  ( yym )  p  end  a  correction  quantity 

<*m  *  (wm)p  +  A"m 


For  the  potential  flow 


1*?^  +  _*?nip  *  0 


(8-103) 
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Further 


(*a»> 

1  dr  ' 


on  dr  dm  '  dr 


dw. 


m 


COS  a  +  ;■ K|  sin 

o  m 


Introducing  all  this  into  equation  (B-101) 


dAw, 


m 


Aw, 


m 


too  £, 


n 


cos  a 


*  {  2"  7* 

%  v  r  9 


-  tan  ft  [^(s.na  * 


sin  a 


dAwm 

dm 


-«•  -4“*  -  ^— >  +  2  77<"m>p 


(B-IOI) 


WA-DC  TR5«-589 


155 

CONFIDENTIAL 


CONFIDENTIAL 


In  this  equation  we  replace 


on 


by 


In  a  mixed  flow  impeller,  it  happen?  that  the  angle  a  at 
any  point  is  very  nearly  equal  to  the  angie  a  ,  taken 
at  the  same  Z  value.  Therefore 


Jilt*  —  a 
‘-os  at 


sinjj. 

sin  a. 


i.O 


tan*  fit 


_rl 

r«* 


)  +  Awm 


p  «<■"*/>» 

r 

2W  — 

tan  fit 

(8-105) 

fr-} 

r. 

COI  a» 

+  tan* 

r* 

tan  a,  jt 

tl'OElrs 
dm  ' 

+  tan*  fit 

Tm 

+  2  cSt?  T?(wm>p  *  *<m' A  ton  °» 


This  equation  is  linear,  except  for  the  presence  of  the  hist  terra. 
The  solution  is  carried  out  by  an  iterative  process.  The  last 
term  is  omitted  and  the  remainder  at  the  equation  solved  for 
each  potential  line.  This  permits  the  evaluation  of  the  term  in 
[d(AWm)]p 
dm 

The  potential  line  solution  is  then  repeated  with  all  terms  in¬ 
cluded,  until  convergence  is  met.  The  solution  of  the  linear 
equation  is  carried  out  as  outlined  under  part  A  above,  except 
for  the  evaluation  of  the  aonstant  of  integration  which  now  has 
the  form 


K  *  - 
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.hub  g  p  (2Tf  .  zt )  dn 

>stir  *J±[±  _ 

^hub  g/>  ( 2*rr  *  Zt  )dn 


shr  ./If.dn 
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APPENDIX  C 

EQUATIONS  FOR  THE  AXI-SYMMETRIC 
FLOW  IN  IMPELLER  PASSAGES  - 
EQUATIONS  FOR  THE  NUMERICAL 
SOLUTION  OF  TEE  DIRECT  PROBLEM 


The  bss*.''  ..^ration  lor  the  numerical  solution  of  the  flow  in  impeller  I 
wan  taken  a»  equation  53  of  ref.  52. 

i  •  .  0  »  d  I  /  .  tan  *3  »  ^  .  3*^ 

37*  ~  T1'  ~  17  +  7F~ 


t  * 

r  p 


\~FZ~  ~‘h+  T  -  ~~z~  +  TT  *  F'l  =  ° 

^  5r  ^0f 

Three  particular  variations  of  this  equation  are  of  interest:  (1)  In 
the  particular  case  of  compressible  flow  in  the  inlet,  C- 1  reduces  to: 


a‘* 

dr* 


a  »  +  a*f  +  j_  a£ 
a  r  az£  P  •  dr  dr 


dr  +  P 


(C-3) 

£  ll  -0 


(2)  In  the  case  of  incompressible  in  vise  id  flow  in  an  impeller  with  radial 
elements  and  no  pre-whirl  (in  which  case  comparison  can  be  rnaA»  with 
equation  B-6): 

tan*£  .  a*t  I  ton*/9  » dt  .  ^  ton£ 

co s  o  dr  r  cos  a  dr  di*  cose 


(C-4) 


157 

CONFIflFNTIM 


tfADC  TO  If- 58® 


CONFIDENTIAL 


If  the  ivi'AVinz  relation*  xre  Introduced,  equation  C-4  reduces  to  B-6 

4r_  df/dr 

cwe  '  dn 

d^f  / d  z 

ton  o»Mr/dz*,»  -3^757  ' 

- 'rJffiW 

,  -[1§7|,*(l71*1% _ . 

'  ,2±m-?  £±  *±  i£.  +  i!±  (Al-t* 

»-*  **  '  drdl  ii  it  Ar*  *  * 


jam 


“M 


id±. td**  +oii  d±  dj±  .  rd»%«d,» 
•di1  dp  _  dr  dz  drd  z  vdr  '  dr*  _ 


(4if  +  (-!■£)« 

1  d  z  dr' 


f) 


a*  at  dt  *  r, 

dT  17  (  dT>~  gHilii] 


I*  +  f  ~  1*  f  J 


Eramn 

PIDH 


)  + 


nff  )*  *  <#?]* 

Equation  C-4  Is  of  interest  in  that  its  solution,  carried  out  siaooltaneoualy 
with  an  exact  expression  for  the  density,  yields  a  solution  directly  comparable 
to  the  quasi-compress  lble  procedure  as  used  In  the  general  design  procedure, 

(3)  The  form  of  equation  C-l,  presently  coded  for  numerical  solution,  makes 
full  allowance  for  compressibility,  but  is  restricted  to  impellers  with  radial 
elements  and  no  pre-whirl.  The  effect  of  losses  Is  Included  to  the  degree  that 
it  affects  the  magnitude  of  the  local  densities;  but  the  effect  of  the  locate 
on  tho  radial  entropy  gradient  is  ejected.  The  clogging  effects  due  to  both 
the  presenca  at  the  blades  and  the  viscous  effects  sure  included. 
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With  these  condition*,  equation  C-l  becomec: 

*  j  +  tqn>£  *  <?If  _L  (t  _  >  djjf 

cos*a  dr*  r  v  ooc*a' dr  dz* 

.  _L  f{  I*  ,  It  5c.  +  £±  «£.]  -2-  =  o  <C' 

*•  cos*  a  dr  dr  dz  dz  **  coso  /Vj 

This  equattoo  (C-5)  must  be  solved  simultaneously  wuh  the  poiyt~opic  density 
relationship 


(C-6) 


For  purposes  of  numerical  computations,  equations  C-5  and  C-6  were  expressed 
on  a  dimensionless  basis,  using  barred  symbols  ( or  dimensionless  quantities 
which  are  defined  as  follows: 


fy 

P 

T 

z 


^/Wc  2  **Pqi 


r  /  h  where  h  is  the  grid  mesh  six*  of  the  flow 

m t. 


2  4 


W 


Wc/2»h%. 
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Further,  the  following  definition  is  made  of  the  over -all  clogging 
factor  «f 


*»  *  *  x  *  Blades 


«  is  viscous  coefllchent  and  is  prescribed  ; 


'Blade® 


I  - 


ton  B  — .  d  o 

We  may  replace  the  ratio  by  its  equivalent  r  jy- 


Writlng  partial  derivatives  with  die  conventional  subscript  notation, 

(W/dz  =  ^z  etc  )  equation  C-5  becomes: 

I  —  P  7  1  y  —  P7  *7  — 

♦Vr  '  T  'V  +  +  *  -T'*T 

-(4?  ♦  +(5Lf+  -=?,F‘  K 


2“  r*  *  «  ? 


tC-7) 


Likewise  In  equation  C-6 


w  *  J5T5 


(C-6) 


T»  MV  A  V- 


wm  *  '/wr*  +  w/ 


-  /TiTTf, 

2»«rh?^ 
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and  where 


cos  0 


\  +  (  r  COS  a 

0  2 


_ *  .  d <t>  .a  (dv/dr)* 

r  (5T’  wtt'  ~ (if/a?)* 


(C-'O) 


Introducing  thin  into  equation  0*0,  there  follows,  alter  simplification: 


■  0+¥< 


k-l  ,  Wc  „  r  .  +1  +  + 

2  'Zw^tf'  'r  “  '  (777?  " 


(C-ll) 


Solution  of  the  flow  in  the  impeller  is  then  reduced  to  the  simultaneous 
solution  of  equation  C-7  and  C-ll. 

It  may  be  recalled  that  from  equation  C-l,  the  term  involving  the  radial 
entropy  gr?  ilent  was  o  mi  tied,  somewhat  arbitrarily,  but  mainly  to  avoid  the 
extra  complicating  features  that  it  introduces.  Thus,  for  a  polytrcpic  process: 

Tds  5  '■  s  -  ft’  RdT 

and  since 


T 

.  ,  +  k-' 

2 

Tor 

Tds 

n  -  k 

(  r«s  - 

dr 

(n- Dkg 

d  w  . 
w  ~ r~  ) 
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An  Idea  of  the  added  complication  may  be  had  by  noting  that  the 
relative  velocity  w  must  be  obtained  from  relations  C-8,  C-9  and 
C-10,  which  themselves  involve  the  density. 
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APPENDIX  D 

EQUATIONS  FOR  THE  CORRECTION  FROM  THE 
TWO-DIMENSIONAL  TO  THE  AXI-SYMMETRIC 
FLOW  UPSTREAM  OF  IMPELLER 


I  Approximate  correction  to  convert  the  velocity  distribution  of  an 
Incompressible  2- dimensional  flow  to  that  of  an  Incompressible  axi- 
symmetrlcal  flow. 

Comparison  is  made  here  of  the  velocity  gradients  along  a 
potential  line  for  a  two-dimensional  flow  passage  and  an  aid- symmetric 
flow  passage  of  identical  meridional  section.  The  standard  correction 
is  taken  from  the  well  known  flows  around  a  sphere  and  a  cylinder,  at 
the  point  of  maximum  velocity  for  a  sphere. 


cm  *  [  jpr  sin  e  +  U  sin  0  \ 


using  standard  classical  notation. 
At  any  angle  9 
dc 


dr 


~  - jr-  U  sin  9 
2  r 


For  a  cylinder 


dc 


m 


o* 


cm  -  U  (  I  +  77  )  sin  9 


-  2  U  -Sj  sin  3 

r3 


And  (^cm/^r)  2-Dim  _  4  r 

(dcm/ar)  3-Dim  3  a 
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Near  the  periphery  - 01 

dr 

(acm/*n)2-Dim  .  _4 

Wcn/«nl3-D""  '  3 


_H?  and  r  5  a 
dn 


We  may  now  extend  the  application  of  this  ratio  to  wd-  tymmetrical 
shapes  for  which  the  meridional  view  Is  other  than  a  circle,  as  for  example 
both  to  the  hub  and  shroud  of  a  compressor  inlet.  And  by  interpolation, 
it  may  be  safely  assumed  that  this  same  ratio  is  applicable  to  Tiny  point 
along  its  potential  line,  between  hub  and  shroud. 

A  check  of  this  principle  was  made  for  a  given  inlet  for  which  both 
the  two-dimensional  and  the  axi- symmetrical  solution  were  available.  The 
results  shown  on  fig.  73  Indicates  the  quantity  4/3  to  be  representative 
of  the  average. 
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BLADE  TO  BLADE  SOLUTION 


Direct  Problem :  Calculation  of  the  blade  loading  fo-  -  rive „ 
geometry.  ~  6*Ten 


-  ^  i  exPre88loin  for  blade  loading  can  be  obtained  by  application 
of  the  laws  of  angular  momentum  to  the  fluid  element  of  figure  74 


Ap  dm  dn  r  =  wcos£  gp  :  A£  dn  d(rcu) 


A £ 
A<£ 


p  w  c 


a  d(fCu) 

B  ~^r 


where  the  barred  quantities  denote  mean  quantities  between  ♦»».»  • 

bl^e  surf  are.  It  is  dear  tnat  these  a«n 

the  corresponding  quantities  of  tha  axi- symmetric  solution. 

Niw  at  a  given  radius  the  total  pressure  is  constant 


Ap  ’  P-,  < wsuct  -  «pr> 

Assuming  that 

=  j  ?<"suct  +  V1 

Ap  =  pw{  w5ucf-  wpr|  =  ?5Aw 
and 

-£r  •  ».?  4^ 

A  dm  (0-1) 
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In  view  of  the  simplifying  asst  rapt! '.ns  Introduced  in  thu  derivation, 
particularly  in  the  handling  of  he  a^itzgc  quantities,  the  resulting 
must  be  regarded  as  an  Irfex  if  th»  olade  loading,  rather  than  a 
tree  velocity  difference.  This  fact  i-o<fs  not  detract  from  its  useful¬ 
ness.  The  relative  velocities  on  tf  »  uressure  and  suction  sides  var>j 
taken  as 


wsuct  s 

.Aw 
"near  2 

"or  * 

w  -  ±-'*- 

"mton 

II  Inverse  problem:  required  bh.^1  ingle  to  produce  a  given  lcad- 


ing 

r  cu  * 

r  (  r  w  - 

wm  ton 

$  > 

d(rcu) 

dm 

*  2  r  aj 

dr  _ 
drn 

dttqnJH 
r"  dm 

-  r  ton 

Q  dWm 

9  dm 

wm  ton  $ 

d) 

dm 


Substituting  Into  equation  D- 1  ir  i  rearranging  terms,  yields 

d(tonfl)  _  ,  2  ru  -  wm  ton  ,}  rj  r 

'  dm 

_  A  w  _  ton  J  dwrp  (J-2) 

Aipcostfrwm  wn  dm 

If  a  tentative  meridional  ve  ocity  iLsvrlbution  has  been  calculated 
and  a  desired  loading  dlstrlbull  n  'Established,  ths  blade  shape  can  Ire 
obtained.  The  equation  is  non-linear  arid  can  best  be  solved  by  yro- 
ceding  in  a  step  by  step  fashion  fn.m  impeller  inlet  to  ejrtt«  For  each 
limit  step,  A,  ton  ft  can  be  computed  since  the  mean  value  of  fi 
Am 

can  be  anticipated. 
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